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ABSTRACT 


A  simple  moment  solution  is  given  for  low  frequency  electromagnetic 
scattering  and  radiation  problems.  The  problem  is  reduced  to  the  correspond 
ing  electrostatic  and  magnetostatic  problems.  Each  static  problem  is  then 
solved  using  the  Method  of  Moments.  The  surface  of  the  perfectly  conducting 
scatterer  is  modeled  by  a  set  of  planar  triangular  patches.  Pulse  expansion 
functions  and  point  matching  testing  are  used  to  compute  the  charge  density 
in  the  electrostatic  problem.  For  the  magnetostatic  current  a  new  set  of 
charge-free  vector  expansion  functions  is  introduced.  The  problem  is  first 
formulated  assuming  the  scatterer  to  be  in  an  unbounded  homogeneous  region. 
Then  the  presence  of  an  infinite  ground  plane  is  incorporated  into  the 
formulation.  Scatterers  of  various  shapes,  such  as  the  circular  disc,  the 
square  plate, the  sphere,  and  the  cube  are  studied.  Special  attention  is 
paid  to  a  conducting  box  with  a  narrow  slot.  The  computed  results  are  the 
scattered  fields,  the  induced  charge  and  current  distributions,  and  the 
induced  electric  and  magnetic  dipole  moments.  These  are  in  close  agreement 
with  whatever  published  data  are  available. 
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Chapter  1 


INTRODUCTION 

Consider  an  arbitrarily  shaped  perfectly  conducting  body  illumi¬ 
nated  by  an  incident  electromagnetic  field  (E1 ,  H1).  The  problem  is  to 
find  the  scattered  field  (ES,  HS)  in  the  Rayleigh  region.  The  Rayleigh  region 
is  defined  [1]  to  be  the  range  of  frequencies  for  which  the  maximum  dimen¬ 
sion  d  of  the  body  is  much  smaller  than  the  wavelength  It  has  been 
shown  by  Rayleigh  [2]  that  this  low  frequency  scattering  problem  can 
approximately  be  solved  by  treating  the  corresponding  electrostatic  and 
magnetostatic  problems  separately.  Also,  it  has  been  established  that  a 
small  scatterer  [1],  (or  a  small  aperture  in  a  perfectly  conducting  screen 
[3])  may  be  approximated  by  radiating  electric  and  magnetic  dipoles.  Exact 
values  for  the  dipole  moments  are  available  in  the  literature  only  for 
simple  shapes,  such  as  the  sphere  and  the  circular  or  elliptical  disc  (or 
aperture)[3]-[6] .  An  extensive  bibliography  on  the  low  frequency  scatter¬ 
ing  problem  is  given  in  [1]. 

Our  work  formulates  the  electrostatic  and  the  magnetostatic  problems 
separately.  An  approximate  solution  is  then  obtained  for  each  problem  using 
the  Method  of  Moments  [7].  In  the  electrostatic  problem  an  integral  equa¬ 
tion  is  solved  for  the  induced  charge  density  on  the  surface  S  of  the  body. 
Then  the  induced  electric  dipole  moment  is  computed  from  the  induced  charge 


distribution.  Similarly,  in  the  magnetostatic  problem  an  integral  equa¬ 
tion  is  solved  for  the  induced  current  density  on  S.  The  induced  magnetic 
dipole  moment  is  computed  using  the  induced  current  density. 


In  principle,  this  scattering  problem  can  be  handled  by  solving 
an  electric  field  integral  equation  for  the  induced  current  on  S.  This 
is  called  an  E-field  solution.  However,  two  static  problems  require,  in 
general,  less  effort.  Furthermore,  although  the  E-field  solutions  give 
accurate  results  in  the  resonance  region,  where  d  is  comparable  to  X,  they 
may  begin  to  lose  accuracy  as  the  frequency  decreases.  Recently,  Mautz 
and  Harrington  [81  nave  discussed  this  problem  in  detail,  and  have  given  a 
new  E-field  solution  for  bodies  of  revolution  which  gives  accurate  results 
in  both  the  resonance  region  and  the  Rayleigh  region. 

The  general  formulation  is  presented  in  Chapter  2.  The  formulation 
of  the  electrostatic  problem  is  known  and  has  been  used  by  various  investi¬ 
gators  [9].  The  approach  is  to  first  model  the  surface  S  by  planar  triangu¬ 
lar  patches  and  then  assume  a  constant  charge  density  on  each  patch.  An 
integral  equation  for  the  charge  density  is  obtained  by  using  the  conditions 
that  the  total  electrostatic  potential  on  S  is  a  constant  and  that  the  total 
charge  is  zero.  This  integral  equation  is  solved  approximately  by  the 
Method  of  Moments. 

The  same  triangulated  model  of  S  is  used  in  the  magnetostatic 
problem.  An  integral  equation  for  the  induced  magnetostatic  current  is 
obtained  by  using  the  condition  that  the  component  of  the  total  magnetic 
field  normal  to  the  surface  S  is  zero  .  This  boundary  condition  is  bor¬ 
rowed  from  the  dynamic  problem,  but  it  holds  for  perfectly  diamagnetic  (y  =0) 
bodies  in  static  fields  as  well.  The  surface  divergence  of  the  induced 
magnetostatic  current  must  be  zero.  This  is  assured  by  introducing  a  set  of 
charge-free  vector  expansion  functions  for  the  current . 

To  compute  the  scattered  fields,  space  is  divided  into  four  regions; 
very  close  to  the  scatterer,  the  electric  field  is  computed  from  the  induced 


charge  only,  and  the  magnetic  field  is  computed  using  the  induced  current 
only.  At  a  distance  large  compared  to  d,  but  small  compared  to  X,  the 
electric  field  is  computed  using  the  static  expression  for  the  field  of 
the  induced  electric  dipole,  and  the  magnetic  field  is  computed  using  the 
static  expression  for  the  field  of  the  induced  magnetic  dipole.  In  both 
of  these  near-zone  cases  the  electric  and  magnetic  fields  are  decoupled. 

For  an  intermediate  region,  where  the  distance  from  the  scatterer  is  compa¬ 
rable  to  X,  the  exact  dynamic  expressions  are  used  to  compute  the  fields 
from  both  dipoles.  In  the  far-zone  region,  where  the  distance  from  the 
scatterer  is  large  compared  to  both  d  and  X,  far-field  expressions  are  used  for 
the  dipole  fields,  and  the  fields  are  coupled. 

In  Chapter  3,  the  general  formulation  described  in  Chapter  2  is  applied 
to  various  scatterers  in  free  space.  Open  and  closed  surfaces  are  considered. 
The  incident  field  is  taken  to  be  either  a  plane  wave  or  is  assumed  to  be 
produced  by  an  oscillating  electric  or  magnetic  dipole  placed  near  the  scat¬ 
terer.  Special  attention  is  paid  to  the  case  where  an  electric  (magnetic) 
dipole  is  placed  normally  (tangentially)  on  a  conducting  sphere  or  on  the 
face  of  a  conducting  cube. 

Chapter  4  formulates  the  scattering  problem  in  the  presence  of  an 
infinite  ground  plane.  In  the  first  part  of  the  chapter,  the  scatterer  is 
assumed  to  be  above  the  ground  plane.  In  the  remaining  part,  the  special 
problem  of  a  conducting  box  on  a  ground  plane  is  considered. 


Chapter  2 


GENERAL  FORMULATION 

In  Appendix  A  we  describe  how  to  approximate  the  problem  of  low 
frequency  electromagnetic  scattering  by  the  corresponding  electrostatic 
and  magnetostatic  problems.  In  this  chapter  we  introduce  the  general 
procedure  for  the  formulation  and  solution  of  these  problems  separately. 

The  main  objective  in  an  electrostatic  problem  is  to  determine  the 
electric  charge  distribution  induced  on  a  perfectly  conducting  iso¬ 
lated  body  when  it  is  placed  in  a  specified  impressed  static  electric 
field.  Once  the  charge  density  is  known,  other  quantities  of  interest, 
such  as  the  induced  electric  dipole  moment  and  the  field  produced  by  the 
induced  charge  ,  can  readily  be  found. 

In  the  magnetostatic  part  of  the  problem  we  are  interested  in  de¬ 
termining  the  current  distribution  induced  on  a  perfectly  diamagnetic 
body  when  it  is  placed  in  a  static  impressed  magnetic  field.  Once  the 
current  density  is  obtained,  other  related  quantities  can  be  solved  for 
easily . 

For  each  problem  we  obtain  a  different  integral  equation  and  use 
the  Method  of  Moments  to  approximate  this  integral  equation  by  a  set  of 
linear  algebraic  equations.  These  algebraic  equations  are  then  solved  by 
the  methods  of  linear  algebra. 

2.1.  The  Electrostatic  Problem 

The  electrostatic  problem  has  been  studied  extensively  and  the  general 
formulation  of  the  problem  in  terms  of  the  Method  of  Moments  is  well-estab¬ 
lished  [7,9].  We  summarize  the  formulation  here  for  the  sake  of  complete- 


ness  and  for  quick  reference.  We  will  refer  to  it  in  the  formulation 
of  the  magnetostatic  problem  in  Section  2.2,  and  also  in  Chapter  4. 
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body . 


Equation  ( 2  —  1 )  is  an  integral  equation  ol'  the  first  kind  for  It  is  to 

be  solved  by  the  Method  of  Moments,  with  the  auxiliary  condition  (2-2). 


2.1.2)  The  Expansion  Functions  for  the  Electrostatic  Charge  Density 

To  solve  for  j,  we  approximate  the  surface  S  by  planar  triangular 
patches  defined  by  a  set  of  nodes  and  edges.  The  vertices  of  triangles  are 
referred  to  as  nodes  and  the  sides  are  referred  to  as  edges  (Figure  2.1). 

We  assume  that  the  charge  density  is  constant  on  each  patch.  In  other  words 
we  let 

N 

-  (r )  =  j  ■:  P  (r)  (2- 

—  ,  n  n  — 


where 


N  is  the  number  of  triangular  patches. 


<7^  ,n  =  1 , 2 , .  .  .  ,N,  are  the  expansion  coefficients  to  be 


determined,  and 


Pn(r)  is  a  pulse  expansion  function  defined  as  follows. 


1 ,  r  in  S 


,n 


P  (r) 

n  ~ 


(2-1 


0,  otherwise 


Here  S  denotes  the  nth  triangle.  With  this  choice  of  expansion  functions. 


a  directly  gives  the  charge  density  on  S  ,  and  (2-2)  reads 
n  n 


Tn  =  0 


(2-1 


n=l 


where  T  denotes  the  area  of  the  nth  triangle.  Substituting  (2-3)  into 
(2-1) ,  we  have 


?  ff  d 

>  a  i  .  -t — 


s  i 

' — ,  r  =  •'*  •  •  V  -  4ttc«j>  (r),  r  on  S 


(2-1 


Fig.  2.1.  A  triangulated  model  of  a  circular  disc.  No.  of  the  triangular 
patches=24,  No.  of  nodes=17,  No.  of  edges=40,  No.  of  interior  nc>des=9. 
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We  could  have  chosen  some  other  types  of  expansion  functions.  For 
example,  we  could  choose  functions  whose  charge  density  is  allowed  to 
vary  on  each  triangle.  This  might  approximate  the  actual  charge  density 
more  closely.  We  could  also  choose  the  expansion  functions  such  that  the 
total  charge  associated  with  each  expansion  function  is  zero.  Then  the 
condition  (2-2)  for  Q  =  0  would  automatically  be  satisfied.  Since  we 
solve  the  electrostatic  problem  as  one  part  of  a  low  frequency  scattering 
problem, we  will  take  Q  =  0. 

The  pulse  expansion  functions  for  the  charge  density  have  been 
successfully  used  by  other  investigators  [9],  and  the  results  are  readily 
available.  We  have  the  additional  advantage  that  for  this  choice  we  can 
make  use  of  the  moment  matrix  obtained  in  the  electrostatic  problem  to 
fill  the  moment  matrix  of  the  magnetostatic  problem  without  additional 
integrations . 

2.1.3)  Testing  Procedure  for  the  Electrostatic  Problem 

For  the  reason  mentioned  above,  and  because  of  its  simplicity,  we 
choose  point  matching,  i.e.,  we  satisfy  (2-6)  at  the  centroid  of  each 
triangle,  obtaining 

N 

I  K  a  -  4tt c  V  =  -  4TTC4)1  ,  m=l,2,...,N  (2-7) 

,  mn  n  o  m 

n=l 
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Here  r_'  is  the  position  vector  of  the  differential  element  of  area  ds*  and 
_rm  is  the  position  vector  to  the  centroid  of  the  mth  triangle  Sm.  The 
integral  in  (2-8)  can  be  evaluated  analytically  and  the  result  is  given 
in  [9]  . 

2.1.4)  The  Moment  Equation  for  the  Electrostatic  Problem 


Using  the  Method  of  Moments,  we  reduce  the  integral  equation  (2-1) 
to  the  set  of  linear  algebraic  equations  given  in  (2-7).  Together  with 
(2-5)  we  can  write  them  in  the  form  of  matrices  as 


K11 

K12  • 

•  *  K1N  ' 

i 

i 

-1  ' 

’ 

°1 

1 

| 

— 4Trect)^ 

K21 

K22 

K2N  ‘ 

f 

-1 

°2 

-4iTe<J)^ 

4 

KN2  • 

•  •  knn  ; 

-1 

°N 

_ _ _ J 

~4iT£<S>n 

T1 

T2  . 

0 

| 

i 

4iTeVj 

Q 

or  as 


' 

[-4ttc({)L] 

m 

[K  ] 
mn 

[-11 

<°o] 

[Tn] 

0 

4tt£V 

o 

Q 

In  matrix  notation  this  is 


(2-10) 


(2-11) 


(2-12) 


where  the  moment  matrix  P  is  the  square  matrix  of  the  order  (N+l)  x  (N+l) 
on  the  left-hand  side  of  (2-11),  "  is  the  column  vector  on  the  left-hand 
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side  of  (2-11),  and  $  is  the  column  vector  on  the  right-hand  side  of 

(2-11).  The  solution  of  (2-12)  gives  both  the  charge  density  and  the 

total  potential  of  the  body. 

All  types  of  electrostatic  problems  can  be  solved  either  by  using 

(2-12)  directly  or  by  modifying  it  slightly.  As  an  example,  if  Q  is 

given  and  is  zero,  then  one  can  obtain  the  charge  distribution,  the 
potential  V  ,  and  the  capacitance  of  the  body  by  using  (2-12)  directly. 

In  fact,  it  is  easy  to  show  that  the  capacitance  of  the  body  is  given  by 

Avr 

C  =  — j-=-“—  (2-13) 

P(N+1,  N+l) 

where  P  ^(N+l,  N+l)  is  the  element  in  the  (N+l)th  row  and  the  (N+l)th 
column  of  the  inverse  of  the  matrix  P.  Another  simple  and  practical  appli¬ 
cation  of  the  electrostatic  problem  formulated  here  is  to  find  the  magnetic 
polarizability  tensor  of  a  small  aperture  of  arbitrary  shape  in  an  infinite 
conducting  screen  [10], [11]. 

As  pointed  out  earlier,  when  we  solve  the  electrostatic  problem  pre¬ 
sented  here  as  a  part  of  a  low  frequency  scattering  problem,  we  take  Q=0 
in  (2-11).  Then,  once  the  charge  distribution  is  solved  from  (2-11),  we 
obtain  the  induced  electric  dipole  moment  according  to  [12] 

Plnd  =  j  o(r')r'  ds'  (2-14) 

S 

Here  _r'  is  the  radius  vector  from  an  arbitrary  origin  to  a  source  point 
on  S.  Since  the  charge  density  is  constant  on  each  triangle,  (2-14)  reads 


N 

> 

n=l 


ff 


ind 


r'  ds' 


(2-15) 


Working  out  the  integral  in  (2-15),  we  have 


pind 


N 


n=l 


n  ' 

n 


,n  n 
r 


(2-16) 


Here  rn  is  the  radius  vector  to  the  centroid  of  the  nth  triangle  and  Tn 
is  the  area  of  the  triangle. 

The  oscillating  electric  dipole  which  is  used  in  computation  of 
the  scattered  field  is  then  given  by 


,  ind  .  „ind 

}_'■  =  JOJ  P 


(2-17) 


where  tu  is  the  angular  frequency  of  the  incident  field.  The  details  of  the 
scattered  field  computations  are  given  in  Section  2.3. 

This  completes  the  general  formulation  of  the  electrostatic  part  of 
the  problem,  assuming  the  body  is  in  an  infinite  homogeneous  medium.  Some 
applications  are  given  in  the  following  chapter.  The  formulation  of  the 
electrostatic  problem  in  the  presence  of  an  infinite  ground  plane  is  given 
in  Chapter  4. 


2.2.  The  Magnetostatic  Problem 

As  pointed  out  earlier,  the  electrostatic  problem  has  been  studied 
extensively  by  using  the  Method  of  Moments,  and  the  results  are  available 
for  an  arbitrarily-shaped  surface.  However,  as  an  approximation  to  the 
low  frequency  scattering  problem,  the  magnetostatic  problem  has  been  applied 
to  only  very  simple  shapes.  Tn  fact,  we  are  not  aware  of  any  such  work 
using  the  Method  of  Moments.  Hence,  in  this  section  we  describe  the  formu¬ 


lation  of  the  problem  in  more  detail. 
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The  magnetostatic  problem  can  be  stated  as  follows.  Given  a  per¬ 
fectly  diamagnetic  body  placed  in  an  impressed  magnetic  field  we 

want  to  find  the  current  distribution  induced  on  the  surface  of  the  body 
such  that  the  normal  component  of  the  total  magnetic  field  is  zero  on 
the  surface  of  the  body  and  that  the  induced  current  density  has  zero 
surface  divergence. 


2.2.1)  The  Integral  Equation  for  the  Magnetostatic  Problem. 

Let  the  total  magnetic  field  be  the  sum  of  the  impressed  field 
and  the  field  H(J)  produced  by  the  induced  surface  current  J.  Then  we 
require  that 


n  *  (Hlmp  +  H(J) )  =  0 


(2-18) 


where  S  denotes  the  surface  of  the  body,  and  n  is  the  unit  (outward) 
normal  vector  to  S.  We  also  require  that 


•  J  =  0 
-s  — 


(2-19) 


where  7^  denotes  the  surface  divergence  operator. 

We  are  here  assuming  a  perfectly  diamagnetic  conductor.  Otherwise 
the  condition  that  the  normal  component  of  the  magnetic  field  is  zero  does 
not  hold  or.  the  conductor  in  a  static  field.  One  can  interpret  J  in  (2-19) 
as  the  zeroth  order  part  of  the  total  current  density  induced  on  an  elec¬ 
trically  small  body  when  the  latter  is  expressed  in  terms  of  a  power 
series  in  k,  where  k  is  the  wave  number  of  the  incident  field  [13,  also 
refer  to  Appendix  A|. 

Equivalently,  we  need  to  solve  for  J  such  that 


H  (.1)  =  -  Himp 
n  n 


on  S 


(2-20) 


and  (2-19)  is  also  satisfied.  The  subscript  n  in  (2-20)  denotes  the  normal 
component.  Equation  (2-20)  represents  an  integral  equation  of  the  first 
kind  for  J  which  will  be  solved  by  the  Method  of  Moments. 

2.2.2)  The  Expansion  Functions  for  the  Magnetostatic  Current 


To  solve  (2-20),  we  approximate  the  surface  S  by  triangular 
patches.  Although  we  are  free  to  use  any  patching  scheme,  we  prefer 
to  use  exactly  the  same  patching  scheme  used  in  the  electrostatic  part. 

This  allows  us  to  use  the  same  integrations  computed  in  the  electro¬ 
static  part  to  fill  the  moment  matrix  for  the  magnetostatic  problem. 

We  let  , 

Nb 

J  =  y  C.J.  (2-21) 

-  . 1-1 
1=1 


where 

Nb  is  the  number  of  the  interior  nodes  for  an  open  surface. 
For  a  closed  surface  Nb  is  equal  to  the  total  number  of  nodes 


minus  one. 


is  the  vector  expansion  function  associated  with  the 
ith  interior  (non-boundary)  node,  and  is  the  ith  expansion 
coefficient  to  be  determined. 

The  expansion  function  associated  with  the  bth  interior  node 
is  illustrated  in  Fig.  2.2.  It  has  the  following  properties: 


i)  On  each  triangle  attached  to  node  L,  the  magnitude  of  is  con¬ 
stant.  This  magnitude  is  equal  to  the  length  of  the  edge  opposite 
node  L  divided  by  the  area  of  the  triangle. 

ii)  On  each  triangle  attached  to  node  L,  is  parallel  to  the  edge 

opposite  to  node  L.  circulates  clockwise  about  the  node  L. 
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In  Fig.  2.2  we  have  five  triangles  attached  to  node  L.  (Note 
that  each  triangle  may  be  in  a  different  plane.)  The  capital  letters 
denote  the  nodes.  On  the  first  surface  S^  attached  to  node  L,  is 
called  J  and  is  given  by 


-LI  T, 


(2-22a) 


Here  FB  denotes  the  vector  drawn  from  node  F  to  node  B,  and  T  is  the 
area  of  the  triangle  S^.  Similarly,  on  the  fifth  triangle  attached 
to  node  L,  J  is  called  J.  and  is  given  by 

Lj  L  J 


j  =  *L 
-L5  TL5 


(2-22b) 


Similar  definitions  hold  for  the  other  triangles. 

iii)  The  component  of  normal  to  an  edge  common  to  two  triangles  is 
continuous.  Hence,  no  line  charge  accumulates  on  any  edge.  This  corn- 
component  of  J  is  given  by 


(JL)  -2/*c 

n 


(2-23) 


where  i  is  the  length  of  the  common  edge  between  two  triangles  attached 
to  node  L. 


iv)  It  is  easy  to  show  that 


•J.  =  0 


(2-24) 


So  that  there  is  no  charge  associated  with  J  .  Hence,  by  virtue  of  (2-21), 
we  see  that  the  necessary  condition  in  (2-19)  is  already  satisfied  by  this 
choice  of  expansion  functions.  Indeed  it  can  be  shown  that  [8, Appendix  B] 


-Lk  ”  -Lk  X  — s'Lk 


(2-25) 
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where  n^  is  th e  unit  normal  vector  to  the  kth  triangle  S^,  attached 
to  node  L.  Here  V  denotes  the  surface  gradient  operator,  and  t, ,  is  a 

—  S  LK 

scalar  function  defined  on  S  ,  .  i,1'  ,  has  a  value  equal  to  2  at  the  node 

L,  and  linearly  decreases  to  zero  at  the  edge  opposite  to  node  L,  ^ 

is  discontinuous,  because  it  suddenly  drops  to  zero  at  the  other  two 

edges  on  S  ,  .  Therefore  7  m  ,  is  impulsive  at  these  edges.  However, 

Lk  — s  Lk 

the  composite  function  ;  ^  consisting  of  ij'^,  con,-^nuous  so 

that  its  surface  gradient  is  not  impulsive.  In  fact,  if  all  triangles 
attached  to  node  L  are  in  the  same  plane,  then  0  defines  a  pyramid-like 
surface  centered  at  node  L,  with  height  equal  to  2.  It  has  as  many  side 
surfaces  as  the  number  of  triangles  attached  to  node  L.  In  Fig.  2.3  we 
have  shown  the  part  of  i;  over  two  triangles,  assuming  both  are  in  the 
xy  plane.  The  encircled  numbers  in  the  figure  denote  the  nodes. 

On  the  first  surface  attached  to  node  L,  i}JL  is  called 
which  is  given  by  the  planar  surface  passing  through  the  points  (x^y^.O), 
(x  ,y  ,0),  and  (0,0,2). 


x(y_  -  yJ  -  y(>s  -  xn)  -  x-,y9  +  ynx 


*L1  =  2 


1  2  7  1  2 


(ylX2  "  Xly2} 


(2-26) 


Similarly,  the  equation  of  the  plane  passing  through  the  points  (x2>y2,0), 
(x  ,y  ,0).and  (0,0,2)  defines  </>l2* 


"\2  '  2 


x(y  .  -  y2>  -  y(x3  -  x2)  -  x2y3  +  y2x 

~~  (y2x3  -  x2Y3) 


(2-27) 


J  =  n  *  V  iJj  ,  =  z  *  V  i]; 

— Ll  — Ll  — syLl  —  's  LI 


(y2  "  yl} 


(x2  -  x^) 


(yl  X2  “  xiy2)/2  ~  (ylX2  "  xiy2)/2 


(2-28) 


^  / 
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where  x,  £,  and  Z  denote  the  unit  vectors  along  the  coordinate  axes.  If 
the  nodes  1,  2,  3,...  encircle  node  L  in  the  clockwise  direction,  as  in 
Fig.  2.3,  then  (y^x,,  -  x^y2)/2  is  the  area  of  the  triangle  S^.  As  a 


result,  (2-28)  reduces  to 


=  12 
— I  i  “  t 

LI 


(2-29a) 


where  12  represents  the  vector  drawn  from  node  1  to  node  2,  and  T  is  the 
area  of  the  triangle  S  Similar  results  hold  for  J  ^ »  A  l3’‘-’ 

However,  if  nodes  1,2,3,...  were  to  encircle  node  L  in  the  counter¬ 
clockwise  direction,  then  the  area  of  the  triangle  would  be  -^y^X2  ~ 
and  (2-28)  would  reduce  to 

3,  .  =  -  (2-29b) 

LI 

It  is  evident  from  (2-29a)  and  (2-29b)  that  J_  always  flows  clockwise  about 
node  L  for  nL1=z.  However, if  then  J_L1,  being  proportional  to  nL1, 

would  flow  counterclockwise  about  node  L.  Therefore,  always  encircles 

in  the  left-hand  sense.  (A  left-handed  screw  whose  axis  is  parallel  to  -LI 
would  advance  in  the  direction  of  when  rotated  in  the  direction  o  f  4r> 

v)  The  expansion  function  can  be  written  as  a  linear  combination  of 
expansion  functions  used  in  f 1 A ]  ,  where  an  expansion  function  is  associated 
with  each  non-boundary  edge.  Here  tT  denotes  the  number  of  triangles  attached 
to  node  L.  Usually  the  number  of  interior  nodes  is  less  than  the  number  of 
interior  edges,  and  hence  we  have  fewer  unknowns  for  the  same  patching 
scheme.  On  the  other  hand,  each  of  our  expansion  functions  is  defined  on 
at  least  three  triangles,  whereas  each  of  those  in  [14]  is  defined  on  a  pair 
of  triangles.  Furthermore,  the  expansion  functions  of  [14]  vary  both  in 
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magnitude  and  in  direction  on  a  triangle,  hut  ours  are  constant  there. 

With  such  a  choice  for  expansion  functions,  on  each  triangle  there 
can  be  up  to  three  expansion  functions  associated  with  the  three 
nodes  of  the  triangle.  On  a  triangle  which  has  a  boundary  edge  only  one 
expansion  function  can  exist.  In  patching  the  surface  S,  care  should  he 
taken  so  that  all  triangles  have  at  least  one  interior  node.  (For  closed 
bodies  this  problem  does  not  exist  and  there  exist  exactly  three  expan¬ 
sion  functions  on  each  triangle,  except  on  those  triangles  attached  to  a 
specific  node  where  there  are  exactly  two  expansion  functions.) 

vi)  It  can  be  shown  that  expansion  functions  defined  on  an  open 
surface  form  a  linearly  independent  set.  For  closed  surfaces  the  number  of 
nodes  Nn  is  equal  to  number  of  interior  nodes.  The  set  of  Nn  expansion 
functions  forms  a  linearly  dependent  set,  and  when  we  remove  one  of  them 
(arbitrarily)  we  are  left  with  a  linearly  independent  set.  For  this 
reason  in  (2-21)  is  taken  to  mean  the  number  of  nodes  minus  one  for  a 
closed  surface. 

vii)  The  constant  nature  of  the  current  flow  with  each  basis  function  within 
a  triangle  may  be  at  first  disconcerting.  Certainly  for  a  large  triangle 
modeling  a  curved  section  of  the  scattering  surface,  one  would  expect  the 
direction  and  the  magnitude  of  the  actual  current  to  vary  within  the  triangle. 
In  other  words,  one  might  ask,  "Is  a  superposition  of  the  basis  functions 
within  a  triangle  capable  of  representing  a  prescribed  current  flowing  in 

an  arbitrary  way  within  the  triangle?"  Unfortunately  the  answer  to  such 
a  question  is  "no",  but  we  remind  ourselves  that  this  is  a  low-frequency 
scattering  problem  where  we  do  not  expect  a  fast  variation  of  the  current 
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except  near  the  boundaries,  or  sharp  edges,  or  on  the  regions  very  close 
to  the  source.  Hence  the  problem  can  be  overcome  by  choosing  a  large 
number  of  small  triangular  patches  to  model  such  regions  of  the  scatterer. 

The  same  procedure  should  be  used  in  solving  the  electrostatic  problem. 

The  expansion  functions  for  the  current  are  capable  of  representing  an 
arbitrary  but  fixed  direction  and  an  arbitrary  magnitude  for  the  current 
within  a  non-boundary  triangle.  Near  an  edge  (of  an  open  surface)  the 
component  of  the  actual  current  normal  to  the  edge  must  be  zero  anyway. 

With  these  basis  functions,  the  current  on  a  boundary  triangle  is  guaran¬ 
teed  to  be  tangential  to  the  edge. 

Having  decided  on  the  expansion  functions,  we  now  take  the  usual 
steps  in  the  Method  of  Moments  technique.  Substituting  (2-21)  into  (2-20) 

we  get  b 

N 

H  (  7  C.  J.)  =  -  Hlmp  on  S  (2-30) 

n  .  “ ,  i—i  n 

1=1 

Using  the  linear  property  of  the  operator  H  ,  we  have 
Nb 

7  C.H  . (r)  =  -  Himp(r)  ,  r  on  S  (2-31) 

.  ■  l  ni  —  n  — 

i=l 

Here,  H  . (r)  denotes  the  normal  component  of  the  magnetic  field  at  a 
ni  — 

point  r  (on  S) ,  produced  by  the  current  basis  function  associated  with 
node  i . 

2.2.3)  Testing  Procedure  and  Moment  Matrix  for  the  Magnetostatic  Problem 

We  assume  a  set  u.;,  j=l,2,...,Nb  of  scalar  testing  functions 

such  that  u.  is  non-zero  onlv  over  the  triangles  attached  to  node  j.  Then 
J 

we  define  the  symmetric  product  of  u  with  any  function  f  to  be  the  integral 
of  their  product  over  the  surface  S,  i.e.. 
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(2-34)  as 
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J  f  f 

Z ( j , i)  =  >  ]  u  (r)  n.  (r)  •  H.(r)ds 

m=l  J  J  im  “  "  _1  “ 


(2-37) 


Here,  1L (r )  denotes  the  magnetic  field  produced  by  the  expansion  function 

J.  at  a  point  r,  n.  (r)  is  the  normal  unit  vector  to  the  m  th  triangle  S. 
-l  -  -.im  -  jm 

which  is  attached  to  node  j,  and  u.  (r)  denotes  the  value  of  u.  at  the 

jm  -  j 

point  r  on  S.  .  Now  using  the  relation 
Jm 


H  =  V  '<  A 


(2-38) 


where  A  is  the  magnetic,  vector  potential,  we  have 


t . 

Z (j , i)  =  l  f  u.  (r)n .  (r)  •  (v  x  A.(r))ds 

m=!  JJ  '  ~  -i  - 


(2-39) 


Here  A^(r)  is  the  vector  potential  produced  by  d_  ,  at  the  point  r. 

Now,  if  we  decide  to  choose  u,  to  be  identically  equal  to  if.  in 

jm  .]"> 

(2-25)  and  make  use  of  Stokes'  theorem,  we  arrive  at  the  following: 


Z(j  ,i)  =  -  ), 


'  (f  (B.  (r 

i  JJ  “ 

O 


)  x  V  p.  )  *  A. (r)ds 

-s  jm  -l  — 


(2-40) 


Here  we  used  the  fact  that  the  value  of  0.  is  zero  on  the  boundary  of  the 

J 

surface  S.  composed  of  t,  triangles  attached  to  node  j.  Then  using  (2-25) 
in  (2-40) ,  we  get 


,J  rr 

Z ( j , i )  =  -  >  I.  (r)  •  A .  ( r ) d s 

m=i  JJ  1 


(2-41) 


Since  J.  is  a  constant  vector  which  denotes  the  expansion  function  .1.  on 

M 


the  mth  triangle  attached  to  node  j,  we  can  take  it  outside  the  integral. 


t . 

j  re 

z  (j  ,  i)  =  -  >  .1 .  •  ;  A  .  (r )ds 

,  — j  m  JJ  — 1  — 

m=1  S 

jm 


We  now  approximate  the  integrand  in  (2-42)  by  its  value  at  the  centroid 


of  the  triangle  S.  to  get 
jm 


t . 


Z(j,i)  =  -  f  (.7.  •  A  (r  .  ))  T;_ 


m=l 


-jm  -i  -njm  jm 


(2 


Here  A. (r  .  )  denotes  the  vector  potential  produced  by  J .  at  the  centroid 
— i  — cjm  r  — i 


r  .  of  the  triangle  S.  ,  and  T.  is  the  area  of  S.  . 
-cjm  °  jm  jm  jm 


Assuming  that  the  node  i  has  t  triangles  attached  to  it,  we  write 


A.(r)  as  follows: 

— l  — 


Ai<D  -  if 


J.  (r') 


ir  -  r 


rr  ds’ 


t . 
i 


r-  y 

4tt 


1-1  J 
S 


r  -  r 


7T  ds' 


iZ 


=  ±  l  j 

hr  L  — i  9 

4  £=1  1 


ds' 


r  -  r' 


(2 


Si? 


The  integral  in  (2-44)  is  the  same  as  the  electric  scalar  potential  at  a 
point  1-,  produced  by  a  constant  surface  charge  density  on  the  triangular 
surface  S^.  This  is  the  same  integral  we  had  in  the  electrostatic  problem 
and  its  analytical  evaluation  is  given  in  [9]. 

Substituting  (2-44)  into  (2-43),  we  get 
t .  t . 


Z(.j,i)  5  - 


4  a 


i 


T.  .1  . 

jm  —  jm 


(  I  i1( 

f.=i 


ds’ 

r  . 
-cjm 


(2- 


Comparing  the  integrals  in  (2-8)  and  (2-45),  we  finally  have 
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•b* 


i 


t .  t 

*«.»  S  L  T  J  •  (  l  JjjKIjm,  if)) 

'  m=l  3  3  ?=1  X 


(2-46) 


i 


i 


5 


R 


N 


B 


Here  K(jm,  i£)  is  the  element  in  the  row  corresponding  to  the  surface  S 
and  in  the  column  corresponding  to  the  surface  of  the  matrix  K  used  in 
the  electrostatic  problem  (Eq.  (2-11)).  Hence  the  moment  matrix  P  of  the 
electrostatic  problem  provides  the  integrals  that  are  required  to  fill  the 
moment  matrix  Z  of  the  magnetostatic  problem. 

Returning  to  (2-36) ,  we  have 

'  ffi . 

jm  -  n 


V(j)  =  -  f 

m=l 


ip.  (r)  Himp(r)ds 


(2-47) 


S. 


If  we  approximate  Himp(r)  bv  Himp(r  .  ),  where  r  .  denotes  the  centroid  of 

n  —  '  n  — cjm  -C]m 


the  triangle  S.  ,  then 


V(j) 


-  f  H 


m=l 


lmP(r  .  ) 
n  — cjm  JJ 

S. 


i  p.  (r)ds 
jm  “ 


(2-48) 


The  integral  in  (2-48)  is  equal  to  the  volume  of  the  pyramid-like  object 
of  height  2.  Hence, 


t . 

.1 


V  ( i )  =  -  4  I  T  ■  H  (r  .  ) 
3  L.  im  n  -cjm 
m=l 


(2-49) 


The  moment  equation  (2-33)  with  the  defining  relations  (2-46),  (2-35), and 
(2-49)  is  solved  to  obtain  the  unknown  expansion  coefficients,  and  the 
induced  current  density  J  is  found  by  (2-21).  Once  the  current  distribution 
is  known,  the  induced  magnetic  dipole  moment  can  be  obtained  [12]  from 


M 


ind  _  JL  [ 

=  2  JJ 


J  (r  )ds 


(2-50) 


If  we  have  N  triangular  surfaces  and  denote  by  J  the  total  current  on 


nth  triangle,  (2-501  reads 


,ind 


r  •  J  (r)ds 


n=L  ^h 


L  V  ff  . 

2  •  ■'  '  j]  r  ds 

n=  1 


(2-51 ) 


Here  S  denotes  the  surface  of  the  nth  triangular  patch.  Working  out  the 
integral  in  (2-51),  we  have. 


M 


i  nd  l  c  ,  n  .  n  ,  n 


V  (r°  x  .i)Tr 


(2-52) 


n=l 


Here  r1'  is  the  position  vector  of  the  centroid  of  Sn,  and  Tn  is  the 


area  of  S  . 

The  oscillating  magnetic  dipole  used  in  computation  of  the  scattered 
field  is  then  given  by 


,,„ind  .  x,ind 
K?.  =  .1  oju  M 


(2-53) 


where  w  is  the  angular  frequency  of  the  incident  field,  and  :i  is  the 
permeability  of  the  medium  surrounding  the  scatterer. 

This  completes  the  formulation  of  the  magnetostatic  problem,  assuming 
the  scatterer  is  in  an  infinite  homogeneous  medium.  The  formulation  of  the 
magnetostatic  problem  in  the  presence  of  an  infinite  screen  is  given  in 
Chapter  4. 

Some  applications  of  the  magnetostatic  problem  formulated  here 
will  appear  in  the  next  chapter.  A  simple  and  practical  application  of 
the  problem  is  to  find  the  electric  polarizability  of  a  small  aperture  in 


an  infinite  conducting  screen. 


2.3.  Computation  of  the  Scattered  Field 

s  s 

To  compute  the  scattered  fields  E  and  if  ,  we  divide  the  space  into 
four  regions. 

i)  Very  close  to  the  scatterer  we  compute  E‘s  directly  from  the 
induced  static  charge  distribution  a  as  follows 
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Here  £  denotes  the  position  vector  of  the  field  point,  _r'  denotes  the 

position  vector  of  the  differential  element  of  area  ds',  Sn  denotes  the 

surface  of  the  nth  triangular  patch,  Tn  is  the  area  of  Sn,  a  is  the 

n 

computed  charge  density  on  Sn,  and  rn  denotes  the  position  vector  to  the 
centroid  of  S  . 

s 

Similarly  we  compute  the  scattered  magnetic  field  H  from  the  in¬ 


duced  magnetostatic  current  .1  as  follows: 
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Here  J  denotes  the  total  induced  current  density  on  the  nth  triangular 
patch  S  . 

ii)  At  a  distance  large  compared  to  the  maximum  linear  dimension  d  of  the 


scatterer,  but  small  compared  to  the  wavelength  X,  we  compute  _E  from  the 
induced  static  electric  dipole  Plnd  and  HS  from  the  induced  static  magnetic 
dipole  M*nd.  The  induced  dipoles  are  placed  at  the  origin,  and  the  scattered 
fields  at  a  point  defined  by  the  radius  vector  _r  are  given  by  the  following: 


-/riind  ind 

3  r(P  •  r)  -  P 


E  (r)  = 


4ttc  r' 


H(r )  = 


3  r(M  •  r)  -  M 


4tt  r 


(2-56) 


(2-57) 


where 


(2-58) 


r  =  — 
—  r 


(2-59) 


iii)  At  a  distance  large  compared  to  d,  but  comparable  with  X,  the 
scattered  fields  are  computed  by  replacing  the  scatterer  by  the  oscillating 
electric  and  magnetic  dipoles  of  (2-17)  and  (2-53).  The  exact  dynamic 
expressions  are  used  to  compute  the  fields  produced  by  them  at  a  point  r 
[refer  to  Appendix  B] . 
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where  r  and  r  are  as  defined  in  (2-58)  and  (2-59). 

iv)  For  the  far  field,  we  use  the  far-field  expressions  for  the  dipole 
fields,  neglecting  the  terms  that  vanish  faster  than  1/r,  i.e., 


s ,  ,  e  '^r  ,  .x  ,,  >  .  ind  , ,  ind- 

E  (r)  =  - k  r  x  f  k  ( E  ■  r )  -  ■ .«;  K  • 

—  -  r  •-  —  — 


(2-62) 


s,  ,  e  1  ,  •  ind  ,  T  ind, 

H  (r)  =  k  _r  Y  (k(iC  *  r)  +  -ou  ly  ! 


(2-63) 


Note  that  while  the  near  field  problem  (cases  (i)  and  ( i i ) )  is 
s  s 

"decoupled,"  that  is,  H"  and  If  are  determined  completely  independently, 
the  field  at  large  distances  from  the  scatterer  is  not  decoupled.  In 
general,  the  electric  dipole  contributes  to  both  the  electric  and  t he 
magnetic  far  field,  and  so  does  the  magnetic  dipole.  Expressions  (2-60) 
and  (2-61)  for  the  fields  in  region  iii  are  also  valid  in  regions  ii 
and  iv.  That  is,  if  kr  •'<  1,  (2-60)  and  (2-61)  reduce  to  (2-56)  and 
(2-57),  and,  if  kr  :■>  1,  (2-60)  and  (2-61)  reduce  to  (2-62)  and  (2-63). 


Chapter  3 


SCATTERING  FROM  SMALL  BODIES  IN  FREE  SPACE 


In  this  chapter  we  use  the  procedure  described  in  the  previous 
chapter  to  compute  the  distribution  of  the  charge  and  the  current  in¬ 
duced  on  the  surface  of  a  small  conducting  body  when  it  is  illuminated 
by  an  electromagnetic  wave.  We  assume  the  body  is  placed  in  free  space 

(c  )  and  the  incident  field  is  either  a  plane  wave  or  is  produced 
o  o 

by  an  oscillating  electric  or  magnetic  dipole  placed  near  the  scatterer. 
Given  the  incident  field  (E1  ,HX)  »  we  find  (f1  such  that 


E*  =  -  7<ji1  (3-1) 

o 

where 

E1  =  lim  E1  (3-2) 

~°  k^O 

and  V  denotes  the  gradient  operator.  We  use  (p1  of  (3-1)  as  the  applied 
potential  in  the  electrostatic  problem  (eq.  (2-1)).  Similarly  we  let 


„imp  .  i 

H  =  lim  H  (3-3) 

k-'-O 

and  use  Hlm^  in  (2-18)  to  solve  the  magnetostatic  problem. 

Tf  the  source  of  the  incident  field  is  an  oscillating  electric  dipole 
I£  placed  close  to  the  scatterer,  then  we  only  need  to  solve  the  cor¬ 
responding  electrostatic  problem.  We  obtain  the  electrostatic  problem 
by  replacing  I-  with  a  quasi-static  electric  dipole  P  such  that 

If  =  jo>  P  (3-4) 


Here  is  the  angular  frequency  of  the  incident  wave  in  the  original 


problem.  We  do  not  solve  any  magnetostatic  problem  because  close  to  an 


oscillating  electric  dipole  the  magnetic  field  vanishes  in  the  limit 
k  -*■  0.  That  is,  Hlm^  as  defined  by  (3-3)  is  zero  in  this  case.  Rv  a 
similar  argument  we  solve  only  the  magnetostatic  problem  if  the  source  is 
an  oscillating  magnetic  dipole  Kf  placed  close  to  the  scatterrr.  We  ob¬ 
tain  the  magnetostatic  problem  by  replacing  Kf  with  a  quasi-static  magnetic 
dipole  M  such  that 

K£.  =  j  .u  M  (3-5) 

and  assuming  the  scatterer  to  be  perfectly  diamagnetic. 


3.1.  Plane  Wave  Scatter ing  from  a  Small  Sphere 

The  exact  results  for  this  problem  are  summarized  in  Appendix  A. 

To  apply  the  method  described  in  the  previous  chapter,  we  take  the  follow¬ 
ing  steps: 

i)  For  convenience  the  sphere  is  assumed  to  be  of  radius  1  meter  and  is 
placed  at  the  center  of  a  spherical  coordinate  system. 


ii)  To  approximate  the  spherical  surface  with  triangular  patches,  the 

surface  is  first  divided  uniformly  into  N,  slices  in  the  ^-direction,  then 

<P 

each  of  these  slices  is  divided  into  +1  patches.  The  poles  (0=0°  and 
0  =  180°)  are  taken  to  be  two  nodes.  Hence,  in  each  <}>-slice  we  have  two 
triangular  and  N  -1  quadrilateral  patches.  The  addition  of  diagonals  to 
all  quadrilaterals  completes  the  triangulation  of  the  spherical  surface. 

The  resulting  model  of  the  spherical  surface  is  shown  in  Figure  3.1.  We 

have  a  total  of  2N  NI  triangular  patches  modeling  the  surface.  The  total 

(■0 

number  of  unknowns  for  the  electrostatic  problem  is  then  2NJ_N  +  1.  The 

0  9 

number  of  nodes  is  N„N\  +  2.  Therefore,  the  number  of  unknowns  in  the 

n  <j) 
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i 

■ 

i 


magnetostatic  problem  is  +  1.  In  this  grid  scheme  the  triangulated 

sphere  is  always  inscribed  inside  the  actual  sphere,  because  the  nodes  are 
taken  to  be  on  the  actual  surface.  To  get  a  partially  inscribed  model  we 
simply  take  the  nodes  to  be  slightly  outside  the  actual  sphere,  obtaining 
a  better  representation  of  the  actual  surface.  Obviously  it  is  impossible 
to  patch  a  spherical  surface  exactly  by  using  a  finite  number  of  triangular 
patches.  However  the  triangular  patches  are  better  than  the  rectangular 
patches  in  terms  of  adequately  modeling  an  arbitrary  surface,  and  they  are 
better  than  quadrilateral  patches  in  terms  of  computer  description. 

For  an  axially  incident  plane  wave  described  by 


we  have 


120'7  e“jk2  x 

(V/m) , 

-  j  kz  x 
e  y 

(A/m) 

-  1  20ii  x 

(V) 

(3-6) 

=  1  V 

(A/m) 

(3-7) 

Figure  3.2  shows  the  computed  current  density  for  this  case.  Also  shown 
in  the  figure  is  the  exact  magnetostatic  result  for  comparison.  The  agree¬ 
ment  is  good.  As  we  noted  earlier,  the  computed  current  is  a  constant  i 

vector  on  each  triangular  patch.  Hence  for  4>  =  0° ,  as  we  change  0,  varies 
on  the  triangular  surface.  However  for  <f  =  90°,  =  -  J^,  which  is  fixed 

on  each  patch.  This  is  the  reason  why  we  have  steps  in  the  curve  for  J  of 
Figure  3.2.  To  observe  the  influence  of  the  patching  scheme  on  the  results, 
we  changed  the  direction  of  F1  and  I]1.  The  change  in  the  results  indicated 


that  the  influence  of  the  triangulation  scheme  is  negligible.  Figure  3.3 
shows  the  computed  charge  density  for 
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6  DEGREES 


Fis-  3.2.  Induced  current  distribution  on  a  sphere  due  to  an  axially 
incident  plane  wave.  No.  of  unknowns  =  17. 
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Qi  =  -  120-ttz  (V)  (3-8) 

The  exact  result  [15]  is  also  shown. 

imp  ^ 

For  an  impressed  magnetic  field  II  =  1  z_  (A/m)  ,  the 

O 

magnitude  of  the  exact  current  at  0  =  90  is  1.5  (A/m),  and  the 
computed  value  is  1.46  (A/m). 

With  an  inscribed  model  of  the  surface,  the  computed  dipole  moments 

are  less  than  the  exact  values.  As  we  slightly  increase  the  radius  of  the 

triangulated  sphere  to  have  a  partially  inscribed  model,  we  notice  that  the 

induced  charge  and  current  distributions  do  not  change,  but  the  induced 

dipole  moments  get  closer  to  the  exact  values.  This  is  to  be  expected,  since 

the  dipole  moments  are  proportional  to  the  cube  of  the  radius,  whereas  the 

current  and  charge  density  do  not  depend  on  the  radius.  Table  3.1  compares 

the  computed  dipole  moments  with  the  exact  values.  Note  that  as  the  radius 

of  the  model  surface  is  increased  from  1.0  to  1.01  (or  1.02),  the  computed 

3  3 

dipole  moments  increase  by  a  factor  (1.01)  (or  (1.02)  ).  Table  3.2  shows 
the  progression  of  the  solution  as  the  number  of  unknowns  is  increased. 

Since  the  scattered  far  field  is  completely  defined  by  the  induced 
dipoles,  we  do  not  need  a  pattern  plot.  From  the  result  we  have  for  the 
induced  dipole  moments  (  the  third  line  of  Table  3.1  )  ,  the  relative  error 
between  the  exact  and  computed  far  field  will  be  at  most  three  percent. 

3.2.  Dipole  Near  a  _Smal 1  Sphere 

Figure  3.4a  represents  an  arbitrarily  oriented  oscillating  dipole 
near  a  conducting  sphere.  The  exact  results,  together  with  some  low  fre¬ 
quency  approximations  for  this  scattering  problem,  are  summarized  in 
Chapter  10  of  [16].  We  first  consider  a  radially  directed  electric  dipole 
near  the  sphere.  Next  we  solve  the  problem  of  a  magnetic  dipole  placed 
tangentially  on  the  sphere. 
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Table  3.1.  Comparison  of  the  computed  and  exact  dipole  moments  for  a 
sphere,  plane  wave  incidence.  N0  =  8,  N ,  =  12,  No.  of  unknowns  in 
magnetostatic  problem  =  97,  No.  of  unknowns  in  e] ectrostat  ic  problem  =  193. 
Exact  area  =  An,  | K1 |  =  120c. 
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'„ind 
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1.1680 xi 0~8  0.49 

Table  3.2.  Progression  of  the  computed  dipole  moments  for  a  sphere,  with 
plane  wave  incidence,  as  the  number  of  unknowns  is  increased.  Radius  of 
patched  sphere  is  taken  to  be  lm. 
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3.2.1)  Electric  Dipole  Near  a  Conducting  Sphere 


Figure  3.4b  shows  a  radially  oriented  oscillating  electric  dipole 
l£  on  a  perfectly  conducting  sphere.  In  the  limit  ka  -►  0,  it  is  stated 
in  [16]  that  the  induced  dipole  moment  is  twice  that  of  the  source  dipole. 
In  fact,  the  first  term  in  the  exact  expression  for  the  total  far  field 
is  the  field  of  3I_^  placed  at  the  origin  (Eq.  (6-118)  of  [4]).  Hence  the 
far  field  can  be  computed  by  replacing  the  system  in  Fig. 3.4b  by  an 
electric  dipole  of  moment  3L#_  placed  at  the  origin.  The  moment  of  this 
dipole  is  independent  of  the  radius  of  the  sphere,  as  long  as  ka  is  very 
small . 

To  get  more  information  about  the  induced  dipole  moment  and  the 
induced  charge  density,  we  used  image  theory  to  solve  the  electrostatic 
problem  depicted  in  Figure  3.5a.  Here  we  have  a  radially  oriented  static 
electric  dipole  near  a  perfectly  conducting  sphere.  To  apply  the  method 
of  images,  we  assume  that  the  dipole  is  composed  of  two  point  charges 
of  opposite  sign  displaced  a  small  distance  d  such  that 


lim  Qd  =  P  (3-9 

d-K) 


The  values  and  the  positions  of  the  image  charges  are  shown  in  Figure  3.5b 
is  used  to  maintain  charge  neutrality.  In  the  limit  d  -*■  0,  we  have  the 
following  results 


2  2  2  2 
Z^Zi  -  5a  )  +  a(Xj  +  3a  )  cos  ') 


2  2  5/2 

(7. .  +  a  -  27  a  cos  D) 


(3-10) 
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1  zi  ’  2  "  (z,  +  d) 

Q3  s  ”  ( Q|  +  Qg  ) 


Fig.  3.5.  a)A  radially  oriented  static  electric  dipole  near  a  conducting 
sphere.  b)An  approximate  representation  of  the  problem  in  (a)  ,  in  terms 
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PLnd  =  2(~)3  P  (3-12) 

Ll 

Here  o  is  the  induced  charge  density  on  the  sphere,  V  is  the  potential 
of  the  sphere,  and  P3nd  is  the  induced  dipole  moment.  Notice  from  (3-12) 
that  for  >•  a  the  induced  dipole  moment  is  twice  that  of  source  dipole, 
independent  of  the  radius  a. 

To  compute  an  approximate  solution  for  the  electrostatic  problem 
of  Fig.  3.5a  we  patch  the  spherical  surface  non-unif ormly .  By  -his  we 
mean  the  density  of  triangular  patches  is  increased  in  the  region  close  to 
the  source.  Since  the  incident  potential  varies  rapidly  in  that  region, 
the  induced  charge  density  is  expected  to  do  so  also.  The  case  where 
the  dipole  is  just  on  the  sphere  (Z^  a)  is  treated  later.  For  Z^  >  a 
we  have  taken 

P  =  4 vzqz  (3-13) 


So  that  the  incident  potential  is  given  by 


z  -  Z1 

((z  -  Zx)2  +  X2  +  y2) 3/2 


(3-14) 


Figure  3.6  shows  the  induced  charge  density  on  the  sphere  when  a  =  2 
meters  and  7.  =  3  meters.  Also  shown  in  the  figure  is  the  exact  result 

of  (3-10)  for  comparison.  The  agreement  is  good.  Figure  3.7  shows  the 
variation  of  the  computed  potential  of  the  sphere  as  the  dipole  moves 
away  from  the  sphere.  The  exact  result  of  (3.11)  is  also  shown.  The 
agreement  between  both  results  is  again  good.  Figure  3.8  shows  the 
variation  of  the  induced  dipole  moment  with  the  location  of  the  source 


dipole.  The  exact  result  of  (3-12)  shows  very  good  agreement  with  the 


computed  one. 


2 


No.  of  unknowns  =  193. 
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Fig.  3.8.  Variation  of  the  induced  dipole  moment  with  the  distance  Z^. 
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No.  of  unknowns  =  193. 
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When  the  dipole  is  just  on  the  sphere  (i.e.,Z^  •+  a),  the  incident 
potential  is  very  large  and  is  changing  rapidly  over  the  region  close  to 
the  dipole.  Hence,  the  direct  use  of  the  point  matching  scheme  does  not 
give  good  results.  To  overcome  this  difficulty  due  to  the  singularity 
in  the  incident  potential,  we  solve  an  equivalent  problem.  Here  another 
dipole  is  placed  just  inside  the  sphere,  below  the  original  source  dipole. 
Figure  3.9  shows  the  configuration.  It  is  clear  that  the  total  potential 
produced  by  the  original  and  the  pseudo-image  dipole  will  be  approximately 
zero  on  a  small  region  defined  by  0q  on  the  spherical  surface  (Fig.  3.9b). 
The  total  incident  potential  outside  this  small  region  will  be  twice  that 
produced  by  the  original  source  dipole,  and  will  vary  slowly.  The  in¬ 
duced  charge  distribution  in  this  equivalent  problem  is  obviously  dif¬ 
ferent  from  the  one  in  the  original  problem.  The  total  potential  at  any 
exterior  point  will  be  the  same  in  each  case,  that  is 


(P)  +  <Ka) 


(3-15) 


or 


<j>T  =  1  ( 2P )  +  <j>(o  )  =  24)1  (P )  +  4>(a  ) 


(3-16) 


Here  <j>  denotes  the  total  potential  at  an  exterior  point,  O  and  a  are  the 

e 

charge  distributions  in  the  original  and  equivalent  problems  respectively, 

♦*(P)  denotes  the  incident  potential  due  to  an  electric  dipole  _P,  and  <J>(o) 
and  <(>(ae)  denote  the  potentials  produced  by  the  induced  charges  in  the  original 
and  equivalent  problems  ,  respectively.  The  induced  dipole  moments  are  related 


as  follows: 
,ind 


pind  =  pind  +  p 


(3-17) 


ind 

is  the  induced  dipole  in  the  original  problem,  P  is  the 

— e 


where  P^ 

induced  dipole  in  the  equivalent  problem  and  P_  is  the  source  dipole  in  the 
original  problem. 


We  assumed  0q  to  be  five  degrees  and  patched  the  dashed  region  of 
Fig.  3.9b  by  twelve  triangular  patches.  Hence,  the  twelve  elements  of 
the  excitation  vector  of  (2-12)  corresponding  to  these  triangles  are 

A 

set  to  zero.  For  a  =  2  meters  and  P  =  Anr  ^z,  the  computed  results  are 
as  follows: 

Plnd  =  2.01125  P 

V  =  -  0.2695 
o 

The  exact  results  are  given  by  the  following 

pdnd  =  2  P 

V  =  -  0.25 
o 

The  relative  errors  in  the  computed  results  are  0.56  and  7.8  percent 
respectively . 

As  pointed  out  earlier,  the  exact  result  for  the  problem  shown  in 
Fig.  3.4b  is  given  in  T4],  Figure  3.10  compares  the  total  computed  and 
exact  fields  at  various  distances.  The  computed  field  is  produced  by  an 
oscillating  electric  dipole  of  moment  3IT  at  the  origin.  The  diameter  of 
the  sphere  is  taken  to  be  2  meters  and  the  wavelength  is  50  meters.  The 
exact  result  is  obtained  by  taking  the  first  seven  terms  of  the  eigen¬ 
function  series  given  in  (6-120)  of  [4].  Note  that  the  fields  in  the 
figure  are  normalized  with  respect  to  the  same  constant.  That  is,  the 
actual  fields  are  |nlt'/4'jk|  times  the  values  shown  in  the  figure.  It 
is  seen  that  even  at  a  distance  equal  to  a  half-wavelength  the  agreement 


between  the  two  results  is  very  good. 
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Figure  3.11  con  .'-ares  the  far  fields  for  two  different  radii.  The 
radius  of  the  sphere  is  taken  to  be  2  meters  or  5  meters.  The  wavelength 
is  fixed  at  50  meters,  and  the  source  dipole  Ij  is  assumed  to  be  on  the 
sphere.  For  both  radii  the  exact  induced  dipole  is  2lj..  Hence  the  far 
field  computed  by  replacing  the  system  of  Fig.  3.4b  with  31 l  placed  at  the 
origin  will  be  the  same  whether  the  radius  is  2  meters  or  5  meters.  This 
is  why  we  have  only  one  curve  for  the  dipole  field  in  the  figure.  For  the 
exact  far  field,  the  first  ten  terms  of  the  eigenfunction  series  are  summed 
for  a  =  5  meters  and  eight  terms  are  taken  for  a  =  2  meters.  The  curves 
shown  in  the  figure  are  normalized,  that  is,  the  true  values  are 

r|T P  e  J  /r  times  the  values  shown  in  the  figure.  It  is  seen  that  if  the 
diameter/wavelength  ratio  is  more  than  0.2,  the  dipole  representation  is 
not  good,  even  for  the  far-field  computations. 

3.2.2)  Magnetic  Dipole  on  a  Small  Sphere 

Figure  3.12a  represents  an  oscillating  magnetic  dipole  K£  tangen¬ 
tially  placed  on  a  small  conducting  sphere.  It  is  stated  in  chapter  10  of 
[16]  that  the  induced  dipole  moment  for  this  problem  (in  the  limit  as  ka  0) 
is  KV/2.  (The  problem  depicted  in  Fig.  3.12a  may  represent  an  equivalent 
problem  for  computation  of  the  field  transmitted  through  a  short  and  narrow 
slot  on  the  conducting  sphere.)  To  solve  the  corresponding  magnetostatic 
problem,  we  use  the  technique  described  earlier,  that  is,  we  place  a 
pseudo-image  magnetic  dipole  just  below  the  original  one, inside  the  sphere, 
and  patch  the  region  near  these  dipoles  densely.  If  Gq  is  small  enough, 
then  the  normal  component  of  the  total  magnetic  field  produced  by  both 
dipoles  is  approximately  zero  on  the  dashed  region  of  Fig.  3.12b.  Hence 
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Fig.  3.11.  The  far  field  for  two  different  radii  in  the  problem  of  Fig.  3.4b. 


the  element  of  the  excitation  vector  ^  in  (2-49)  corresponding  to  the 
node  at  the  pole  (9  =  0°)  is  set  to  zero,  and  there  is  no  contribution 
from  the  triangles  in  the  dashed  region  to  the  elements  of  (2-49)  corre¬ 
sponding  to  the  nodes  on  the  boundary  of  the  dashed  region.  We  failed 
to  obtain  reasonable  answers  even  to  this  equivalent  problem.  We  believe 
this  is  due  to  the  fact  that  the  singularity  of  the  magnetic  field  of 
the  magnetic  dipole  is  of  the  third  order.  Hence  the  variation  of  the 
incident  magnetic  field  is  too  fast  and  the  approximation  of  (2-47)  by 
(2-48)  is  not  valid  in  this  case.  (In  the  electric  dipole  case  the 
singularity  of  the  incident  potential  was  of  second  order,  and  we  could 
obtain  reasonable  results  for  the  equivalent  problem  shown  in  Fig.  3.9b.) 

Notice  that  the  singularity  of  the  electric  field  produced  by  a  magnetic 
dipole  is  of  the  second  order.  Hence  we  use  Maxwell's  equations  and  some 
vector  calculus  to  express  the  excitation  vector  for  the  magnetostatic 
problem  (Eq.(2-47))  in  terms  of  the  incident  electric  field.  From  (2-47), 
we  have 


V(j) 


f  <r)  n  •  H  ”p(r)  ds 

m=l  Jl  " Jm 


V 


xEimp 

jwp 


ds 


(3-18) 


For  the  dipole  of  Fig.  3.12a,  we  have 
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Eimp  = 
where 


and 


Me-jkr  (j|+i)(i(z-a)  -£  y) 
r  r 


2  2  2 
r  =  x  +  y  +  (z-a) 

jcayM  =  K£ 


(3-19) 


(3-20) 

(3-21) 


Substituting  (3-19)  into  (3-18)  and  taking  the  limit  k ->■  0,  we  have 

t . 


v«>  ^  t 

m=l 


T. 


jm  —  jm 


y  (z  .  -  a)  -  z  y  . 

(  -  C3JD _ ~  cjm  n 

V  3 

r  . 
cjm 


(3-22) 


Here  z  .  and  y  .  represent  the  coordinates  of  the  centroid  of  the  tri- 
cjm  cjm 

angle  S.  ,  and  r  .  is  the  distance  from  the  dipole  to  this  centroid. 
e  jm  cjm 

Using  (3-22)  instead  of  (2-49)  for  the  excitation  vector,  we  solve  the 

magnetostatic  problem  shown  in  Fig.  3.12b.  With  dipole  moment 

M  =  0.4ttx,  a=l  meter,  and  6  =5  degrees  (Fig.  3.12b)  the  computed  ratio 

—  —  o  ' 

of  the  induced  to  original  dipole  moment  is  0.4889.  The  exact  result  for 
this  ratio  is  0.5.  The  relative  error  is  2.2  percent.  The  number  of  un¬ 
knowns  used  for  this  problem  is  97.  We  have  also  observed  that  the  computed 
ratio  does  not  depend  on  the  radius  of  the  sphere. 


3.3.  The  Circular  Disc 

Figure  3.13a  represents  a  conducting  circular  disc  lying  in  the  xy 
plane.  The  other  parts  of  the  figure  show  various  types  of  electromagnetic 
illumination  of  the  disc.  The  disc  problem  has  been  studied  extensively 
and  a  summary  of  the  results  can  be  found  in  Chapter  14  of  [16]. 

Figure  3.14  shows  the  induced  magnetostatic  current  density  on  the 
disc  due  to  an  incident  plane  wave,  as  shown  in  Fig.  3.13b.  The  radius  of 


the  xy  plane.  d)A  magnetic  dipole  perpendicular  to  the  plane  of  the  disc. 


EXACT 

COMPUTED 
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the  disc  is  taken  to  be  1  meter,  and  the  number  of  unknowns  used  to 
compute  the  current  density  is  81.  The  exact  result  is  taken  from  [17]. 

It  is  seen  that  the  computed  result  is  a  good  step  approximation  to  the 
exact  one. 

Tables  3.3  to  3.5  summarize  the  results  concerning  the  induced 
dipole  moments  for  the  disc  problem.  The  exact  results  are  obtained 
using  the  equations  (14.185) , (14.187) , (14. 188) , and  (14.193)  of  [16]. 

3.4.  The  Square  Plate  and  the  Bent  Plate 

Figure  3.15  shows  the  induced  current  density  on  a  square  plate  of 
1  meter  side  length,  for  two  different  incident  plane  waves.  The  plate  is 
placed  in  the  xy  plane.  Also  shown  in  the  figure  is  the  result  obtained  using 
the  E-field  solution  program  written  by  Rao  [14],  where  we  let  the  wave¬ 
length  be  56  meters.  It  is  seen  from  the  figure  that  the  magneto¬ 
static  part  is  dominant  in  the  total  induced  current,  even  for  0  =  20 
degrees.  At  0  =  0  degrees  (normal  incidence)  the  magnetostatic  current 
vanishes. 

Figure  3.16  shows  the  induced  magnetostatic  current  density  on  a 
bent  square  plate  of  side  length  1.5  meters.  The  bend  is  located  at  a 
distance  of  one  third  the  width  from  an  edge,  and  a  plane  wave  with  the 
magnetic  field  polarized  perpendicular  to  the  larger  section  of  plate  is 
incident.  The  smaller  plate  section  is  bent  through  an  angle  of  50 
degrees  away  from  the  direction  of  propagation,  which  is  parallel  to  the 
larger  section  of  the  plate. 


Table  3.3.  Induced  dipole  moments  for  the  problem  of  Fig.  3.13b 

E*  =  377  (V/m) ,  H*  =  1  (A/m);  N,  N^  =  No.  of  unknowns,  exact  M^11 
x  z  — 

-2.66  z  .  exact  p*nc*  =  0.178  x  10  ^  x. 


Computed  M  Relative  %  Error 


2.2853 

2.3407 

2.5162 


Computed  P  Relative  %  Error 


0.1539  *  10 


0.1574  x 10 


0.1685  x  10 


Table  3.4.  Induced  electric  dipole  moment  for  the  problem  of  Fig.  3.13c 
No.  of  unknowns  =  109,  source  dipole  £  =  4ire  0  at  (x  ,y  ,z  ). 


Computed  P*“~ 


Exact  P 


%  Error 


0.335 

10 

0.3544  x  10 

0.602 

io-n 

0.63404  x  10 

Table  3.5.  Induced  magnetic  dipole  moment  for  the  problem  of  Fig.  3.13d 
No.  of  unknowns  =  49,  source  dipole  M  =  4tt£  at  (x_ ,  y_,  z_) 


Computed  M 

2 

0.671 

-0.456 


Exact  M 

z 

0.678 

-0.466 


/  Error 


'  1 1 . 


Fig.  3.15.  Induced  current  on  a  square  plate  due  to  two  different 

incident  plane  waves. 


Fig.  3.16.  Distribution  of  the  magnetostatic  current  on  a  bent  square  plate 


No.  of  unkowns  =  77. 
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3.5.  The  Small  Conducting  Cube 


3.5.1)  Plane  Wave  Incidence 

Figure  3.17a  represents  a  conducting  cube  of  side  length  O.lA 
illuminated  by  a  plane  wave.  Figure  3.18  shows  the  induced  current 
density  on  the  right  face,  computed  using  the  magnetostatic  formulation. 
Figure  3.19  shows  the  induced  current  on  the  top  face.  Also  shown  are 
the  results  obtained  by  using  Rao's  [14]  E-field  solution.  Close  agree¬ 
ment  between  the  two  results  is  not  expected,  since  the  maximum  dimension 
of  the  cube  is  0.17A  in  this  case.  Our  result  also  has  some  similarity 
with  the  result  of  [18]  obtained  using  rectangular  patches,  but  a  pre¬ 
cise  comparison  was  not  possible.  Although  the  current  distribution  is 
not  in  close  agreement  with  the  two  results  mentioned  above,  the  computed 
radar  cross  section  does  agree  very  well  with  published  data.  If  of 

Fig.  3.17a  is  taken  to  be  120 7T  (V/m),  then  the  computed  induced  dipoles 

i.nd  ~~~J  /\ 

for  the  cube  of  side  length  one  meter  are  given  by  P  =  0.115x10  x  (C-m  ) 

ind  A  O 

and  M  =  -1.57  z  (A-m).  The  radar  cross  section  obtained  using  these 

2 

dipoles  as  the  source  of  the  scattered  field  is  0.0031X  .  This  result  agrees 

with  the  experimental  and  computational  results  mentioned  in  [ 18 ] — [ 19 ] . 

In  fact,  we  have  observed  that  the  radar  cross  section  obtained  using  these 
dipoles  agrees  (within  less  than  ten  percent)  with  the  published  results 
[19]  even  for  a  side  length  of  0.2A.  Note  that  in  this  case  the  maximum 
dimension  of  the  cube  is  as  large  as  0.35X. 

Finally,  we  have  observed  that  the  induced  dipole  moments  are  pro¬ 
portional  to  the  volume  of  the  cube.  This  is  expected  ,  since  the  induced 
charge  and  current  distributions  do  not  depend  on  the  size  of  the  cube. 


Fig.  3.17.  a)Plane  wave  normally  incident  on  a  small  conducting  cube.  b)An 
oscillating  magnetic  dipole  on  the  top  surface  of  the  cube,  c)  Equivalent 
magnetostatic  problem  for  part  (b)  ;  A  pseudo-image  magnetic  dipole  is 


placed  just  inside  the  cube,  below  the  source  dipole  M. 


Fig.  3.19.  The  induced  current  on  the  top  face  of  the  cube  of  Fig.  3.17a. 
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We  have  also  observed  that  the  induced  dipole  moments  for  the  cube 
are  slightly  larger  than  the  ones  obtained  for  the  sphere,  pro¬ 
vided  the  incident  field  is  the  same  and  the  volumes  of  the  sphere  and 
the  cube  are  equal. 

3.5.2)  A  Magnetic  Dipole  on  the  Cube 

Consider  a  cube  with  some  internal  sources  and  a  small  narrow  slot 
on  one  face.  The  problem  is  to  compute  the  fields  transmitted  through 
the  slot.  If  the  cube  is  small,  then  an  approximate  equivalent  problem 
is  shown  in  Fig.  3.17b.  Here  we  have  a  magnetic  dipole  KJi  placed  tangen¬ 
tially  on  the  cube  at  the  position  of  the  slot,  and  the  slot  is  shorted. 

To  solve  the  corresponding  magnetostatic  problem,  we  use  the  pseudo-image 
technique  to  obtain  the  problem  depicted  in  Fig.  3.17c.  An  image  dipole 
is  placed  just  below  the  original  one,  inside  the  cube.  Hence  the  normal 
(tangential)  component  of  the  total  magnetic  (electric)  field  produced  by  the 
two  dipoles  is  zero  on  the  face.  At  any  other  point  on  the  cube  surface 
the  field  is  twice  that  produced  by  the  original  dipole  K2  alone.  Because 
of  the  singularity  problem  mentioned  in  Section  3.2.2,  we  use  (3-18)  for 
the  excitation  vector. 

The  center  of  the  cube  is  placed  at  the  origin  of  a  Cartesian  coordi¬ 
nate  system.  The  side  length  is  assumed  to  be  1  meter.  The  point  (x  =  0, 

y  =  0,  z  =  0.5)  defines  the  center  of  the  top  face.  The  source  dipole  is 
^  2 

taken  to  be  M  =  lx  (A-m  ).  Table  3.6  shows  the  computed  induced  dipole 

moment,  as  the  position  (x  ,y  ,z  )  of  the  source  dipole  M  is  changed.  From 

o  o  o  — 

the  table  we  conclude  that  a  narrow  slot  at  the  center  of  a  face  of  the 
cube  transmits  energy  outside  less  efficiently  than  a  slot  closer  to  the 
edges  of  the  face. 
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Finally  we  have  the  following  observations:  i)  As  the  size  of  the 
cube  is  changed  keeping  the  source  dipole  at  the  same  relative  position  , 
the  induced  dipole  moment  does  not  change.  For  plane  wave  incidence  the 
induced  dipole  moments  are  proportional  to  the  volume  of  the  cube.  ii)  When 
a  magnetic  dipole  M  is  placed  tangentially  at  the  center  of  a  face  of  the 
cube,  the  computed  induced  dipole  moment  is  0.15  M.  This  is  about  three 
times  smaller  than  the  induced  magnetic  dipole  moment  if  M  were  placed 
tangentially  on  a  sphere.  To  see  the  effect  of  an  electric  dipole  £  placed 
normally  at  the  center  of  a  face  of  the  cube,  we  compute  the  induced  dipole 
moment.  The  computed  induced  electric  dipole  moment  is  0.65  P.  This  also 
is  about  three  times  smaller  than  the  induced  electric  dipole  moment  when 
is  placed  normally  on  a  sphere. 
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Table  3.6.  Induced  magnetic  dipole  moment  due  to  a  tangential  magnetic 

/v  2 

dipole  M  =  1  x  (A-m  )  on  the  top  face  of  a  cube.  The  top  face  is  at 

z  =  0.5.  M  is  at  (x  ,  y  ,  0.5). 

—  o  o 


X 

o 

r~ 

yo 

Induced  Dipole  Min<^ 

0 

0 

0.15  x 

0.2 

0 

0.18  x  -  0.237  z 

0.3 

0 

0.22  x  -  0.406  z 

0.45 

0 

0.36  x  -  0.716  z 

0 

0.2 

0.164  x 

0 

0.3 

0.204  x 

0 

0.45 

0.507  x 

0.3 

0.3 

0.244  x  -  0.33  £  -  0.382  z 

Chapter  4 


LOW  FREQUENCY  SCATTERING  IN  THE  PRESENCE  OF  A  GROUND  PLANE 

Figure  4. la  represents  an  electrically  small  conducting  body  Bl, 
above  an  infinite  ground  plane,  illuminated  by  an  electromagnetic  source 

51.  Using  the  method  of  images,  we  obtain  the  equivalent  problem  shown 
in  Fig.  4.1b.  Here  we  have  two  bodies  Bl  and  B2  and  two  sources  SI  and 

52.  The  ground  plane  is  removed  and  the  space  surrounding  the  bodies  is 
characterized  by  (e,p).  B2  is  the  mirror  image  of  the  true  body  Bl,  and 
S2  is  the  image  of  the  true  source  SI.  To  use  our  low-frequency  approxi¬ 
mation  technique,  we  assume  that  the  true  body  Bl  is  close  to  the  ground 
plane  in  Fig.  4.1a.  More  precisely,  it  is  assumed  that  we  can  inscribe  the 
true  body  Bl  and  its  image  B2  (Fig.  4. lb)  in  a  sphere  of  diameter  d.  Then, 
provided  d  is  much  less  than  the  wavelength  of  the  incident  wave,  we  can 
approximate  the  low-frequency  problem  shown  in  Fig.  4.1b  by  the  correspond¬ 
ing  electrostatic  and  magnetostatic  problems.  The  method  described  in 
Chapter  2  can  directly  be  applied  to  the  two  body  problem  of  Fig.  4.1b  to 
compute  the  induced  charge,  current  and  dipole  moments.  However  this  leads 
to  unnecessarily  large  matrices.  In  this  chapter  we  describe  how  to  avoid 
these  large  matrices  and  formulate  the  problem  in  terms  of  moment  matrices 
which  have  the  same  sizes  as  if  the  true  body  Bl  were  in  free  space.  The 
elements  of  the  matrices  are,  of  course,  more  complicated. 

4.1.  The  Electrostatic  Problem  in  the  Presence  of  a  Ground  Plane 

The  problem  here  consists  of  finding  the  induced  static  charge  dis¬ 
tribution  on  the  conducting  bodies  Bl  and  B2  of  Fig.  4.1b  when  the  incident 


electric  field  due  to  SI  and  S2  is  replaced  by  its  electrostatic  limit. 


Assume  for  the  moment  that  the  true  body  B1  is  strictly  above  the  ground 
plane  in  the  original  problem,  is  isolated, and  has  an  initial  charge  Q  on  it. 
Then  the  net  charge  on  B2  of  Fig.  4.1b  is  equal  to  -Q.  If  the  total  potential 
produced  by  the  sources  SI,  S2, and  the  induced  charges  and  is  VI  on 
Bl,  then  it  is  -VI  on  B2.  (Here  and  are>  respectively,  the  induced 
charge  densities  on  Bl  and  B2.)  In  other  words  we  have  the  following  con¬ 
ditions  to  be  satisfied  by  the  induced  charges: 


.Bl 


<J)S1(r)  +  (J>S2  (r)  + 


■  o1(r’)  ds' 

•  a2(r^')ds' 

4tt£  _r-r '  |  + 

4tt£  r-r* 

(4- 


Bl 


B2 


(“VI,  r  on  S 


B2 


ai(r)ds'  =  - 


a2(r)ds'  =  Q 


(4- 


B1 


B2 
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Here  <j>  (r)  and  <j>  (r)  denote  the  potentials  produced  at  a  point  £  by  the 
true  source  SI  and  the  image  source  S2,  respectively.  The  integrals  in  (4-1) 
represent  the  potentials  produced  by  the  induced  sources,  and  Sg^  and  Sg2 

denote  the  surfaces  of  Bl  and  B2,  respectively.  As  in  Chapter  2,  we 
assume  that  SD1  is  modeled  by  N  triangular  patches  and  the  charge  density 
on  the  nth  patch  is  a  constant  o  .  Then  Sg2,  being  the  mirror  image  of  S^, 
is  represented  by  N  triangular  patches,  and  the  charge  density  on  the  n'th 
patch  (which  is  the  image  of  the  nth  patch  on  S^)  is  Hence,  following 

the  same  steps  used  in  Section  2.1,  we  reduce  the  integral  equation  (4-1) 
to  the  following  set  of  linear  equations: 
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Here  each  subscripted  K  denotes  the  potential  produced  by  a 
constant  charge  density  of  4ne  coulomb  per  square  meter  lying  on  a  tri¬ 
angular  patch  designated  by  the  second  subscript,  at  the  centroid  of  the 
triangular  patch  designated  by  first  subscript.  For  example  K,.^,  denotes 
the  potential  at  the  centroid  of  the  fifth  triangle  produced  by  the  con¬ 


stant  charge  density  over  the7’th  triangle.  The  7  ’  th  triangle  is  on  S  _ ,  and 

dZ 

SI 

is  the  image  of  the  7th  triangle  on  S  .  (j>  denotes  the  potential  pro- 

ol  tn 

duced  by  the  true  source  SI  at  the  centroid  of  the  mth  triangle.  Similarly, 

S2 

4)  is  the  potential  due  to  image  source  S2  at  the  centroid  of  the  mth 
m 


triangle.  Observe  that 


K  =  K  ,  , 
ran  m  n 


(4-4) 


K  ,  *  K  , 
mn  m  n 


(4-5) 


fsl  +  tfs2 

m  m 


..SI  ,S2. 

-  (4>  .  +  <p  . ) 

m  m 


(4-6) 


Then  it  is  obvious  that  (4-3)  represents  the  same  N  equations  written 
twice.  We  approximate  (4-2)  by 


12  N 

T  ai  +  T  02  +  . . .  +  T  aN  -  Q 


(4-7) 


and  combine  (4-7)  with  the  first  N  equations  of  (4-3)  to  obtain  the 
following : 


K11  "  Kl’l  K12  "  K1 ’ 2  ‘ ‘ ‘  K1N 
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In  matrix  notation  this  is 


[PC]  [ZC]  =  [&1  (4-9) 

where  the  moment  matrix  PG  is  the  square  matrix  of  order  (N+l)  x  (N+l)  on 
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the  left-hand  side  of  (4-8),  Fg  is  the  column  vector  on  the  left-hand  side 
of  (4-8) ,  and  $G  is  the  column  vector  on  the  right-hand  side  of  (4-8) .  The 
solution  of  (4-9)  gives  both  the  induced  charge  density  a  and  the  potential 
of  the  true  body  B1  above  the  ground  plane  (Fig.  4- la).  Notice  that  the 
entry  in  the  last  row  and  the  last  column  of  the  inverse  of  the  matrix  PG  is 
equal  to  4+e  times  the  reciprocal  of  the  capacitance  of  the  true  body  above 
the  ground  plane. 

Observe  that  the  number  of  equations  to  be  solved  in  (4-9)  is  the  same 
as  the  number  of  equations  in  (2-12)  where  we  had  no  ground  plane.  However, 
the  matrix  elements  and  the  excitation  vector  elements  in  (4-9)  are  more 
complicated.  Once  the  charge  density  is  computed,  we  use  (2-16)  to  compute 
the  induced  electric  dipole  moment.  Note  that  in  this  case  the  induced 
dipole  radiates  in  the  presence  of  the  ground  plane. 

Figure  4.2a  shows  the  computed  capacitance  of  a  sphere  of  lm  radius, 
the  center  of  which  is  a  distance  d  above  the  ground  plane.  Also  shown  are 
the  exact  results  [12, Eq. (4-42) ] .  The  agreement  is  good.  We  see  that  the  computed 
results  are  less  than  (at  most  by  four  percent)  the  exact  values.  The  reason 
for  this  is  that  the  patched  surface  is  taken  to  be  an  inscribed  model  of  the 
actual  spherical  surface.  As  the  number  of  patches  is  increased,  the  tri¬ 
angulated  model  represents  the  actual  surface  more  closely,  and  the  computed 
capacitance  approaches  the  exact  value.  This  can  be  observed  from  Figure  4.2b, 
which  shows  the  computed  capacitance  versus  the  number  of  patches  for  the 
case  d  =  2  meters. 

We  note  in  passing  that  the  condition  that  the  true  body  B1  be  close 
to  the  ground  plane  is  necessary  onlv  i f  we  use  the  present  procedure  to 
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approximate  a  low  frequency  scattering  problem.  Obviously,  if  we  are 
interested  in  strictly  electrostatic  results,  that  condition  is  no  longer 
necessary. 


4.2.  The  Magnetostatic  Problem  in  the  Presence  of  a  Ground  Plane 

Here  the  problem  consists  of  computing  the  induced  current  distri¬ 
bution  J  on  the  true  body  B1  of  Figure  4.1a  when  the  incident  magnetic  field 
due  to  source  SI  is  replaced  by  its  magnetostatic  limit.  We  again  assume  that 
B1  is  strictly  above  the  ground  plane.  We  first  formulate  the  problem  as  a 
two  body  magnetostatic  problem,  using  the  procedure  of  Section  2.2,  and  then 
use  the  information  provided  by  the  method  of  images  to  reduce  the  number 
of  equations  by  a  factor  of  two.  For  the  two  body  problem  shown  in  Fig. 

4.1b,  the  boundary  condition  to  be  satisfied  is  that  the  normal  component 
of  total  magnetic  field  produced  by  the  sources  SI  and  S2  plus  the  induced 
currents  on  El  and  B2  must  be  zero  on  both  Sn1  and  S„„.  The  condition  that 
the  induced  currents  on  both  bodies  are  charge-free  is  automatically  satis¬ 
fied  by  the  choice  of  expansion  functions,  described  in  Section  2.2.2. 

Assume  that  S  is  patched  so  that  it  has  N*3  interior  nodes.  Then 
d1 

SDO,  being  the  mirror  image  of  S  ,  has  N*3  interior  nodes.  (For  closed 

¥>Z  di 

bodies,  represents  the  total  number  of  nodes  minus  one.)  We  have 
expansion  functions  on  each  of  Sgl  and  S^.  Then,  corresponding  to  (2-33), 
we  have  the  following  set  of  linear  equations  for  this  two  body  problem: 

[  [ZBB]  [ZBI]'  ’  [  TB]  I 

=  [VB]  (4-10a) 

[H]_ 

nziBi  [ziii;  ‘nm 


=  [tfi] 


(4-10b) 
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Here  ZBB,  ZBI,  ZIB,  and  ZIT  are  N  *  N  square  matrices;  IB,  II,  VB,  and 
VI  are  Nb  *  1  column  vectors  .  ZBB(j,i)  is  the  normal  component  of  the  mag¬ 
netic  field  produced  by  the  ith  expansion  function  attached  to  the  ith  node  on 
the  true  body,  tested  over  the  domain  of  the  jth  expansion  function  on  the  true 
body,  as  given  by  (2-34)  and  (2-46).  Similarly,  ZBI(j,i)  represents  the  field 
produced  by  the  image  of  the  ith  expansion  function  tested  over  the  domain  of 
the  jth  expansion  function.  ZIB(j,i)  is  the  field  produced  by  the  ith  expansion 
function  tested  over  the  domain  of  the  image  of  the  jth  expansion  function,  and 
finally  ZII(j,i)  represents  the  field  of  the  image  of  the  ith  expansion  function 
tested  over  the  domain  of  the  image  of  the  jth  expansion  function.  The  ith 
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element  of  IB  is  the  expansion  coefficient  of  the  ith  expansion  function, 

and  the  ith  element  of  II  is  the  coefficient  of  the  image  of  the  ith  expansion 

->■ 

function.  The  kth  element  of  VB  is  the  negative  of  the  normal  component  of 
the  total  impressed  magnetic  field  (produced  by  SI  and  S2)  tested  over  the 
domain  of  the  kth  expansion  function,  as  given  by  (2-36).  Similarly  ,  the 
kth  element  of  VI  is  the  negative  of  the  impressed  field  tested  over  the 
domain  of  the  image  of  the  kth  expansion  function.  Since  Sg2  and  S2  are  the 
images  of  Sg^  and  SI,  and  due  to  (2-46)  and  (2-49),  it  is  not  difficult 
to  see  that 


IB(i)  =  II(i)  i  =  1,2 . Nb, 

VB(k)  =  VI (k)  k  =  1,2,. ..,Nb, 


(4-11) 

(4-12) 


and 

ZBB(j,i)  +  ZBI(j.i)  =  ZIB(j,i)  +  ZII(j,i)  (4-13) 

Hence,  (4-10a)  and  (4-10b)  represent  the  same  set  of  equations.  (This  is 
to  be  expected,  since  satisfying  the  boundary  condition  mentioned  earlier 


on  S  , (4-10a) ,  automatically  satisfies  the  same  condition  on  Sg2  ,  (4-10b)). 
B1 

We  rewrite  (A-lOa)  in  matrix  notation  as 


[ZG]  [I$1  =  [VG] 


(4-14) 


where  ZG  is  a  square  matrix  of  order  N*3  x  such  that 


ZG(j,i)  =  ZBB(j,i)  +  ZBI  (j  ,i) 


(4-15) 


IG  and  VG  are  column  vectors  which  are  the  same  as  IB  and  VB  of  (4-10a) , 
respectively.  The  term  ZBB(j,i)  in  (4-15)  is  equal  to  Z(j,i)  of  (2-46), 
and  ZBI(j,i)  is  given  by  the  following: 


ZBI(j.i) 


t . 


T.  J. 

jm  -jm 


t . 

l 

<  I 

£=1 


K(jm, 


i’ll)) 


(4-16) 


Here  represents  the  value  of  the  image  of  ith  expansion  function  on  the 
5-th  triangle  attached  to  node  i,  and  K(jro,i’£.)  represents  an  element  of 
the  matrix  on  the  left-hand  side  of  (4-3).  The  row  and  column  of  this 
element  are  defined  by  two  triangular  surfaces.  The  row  is  defined  by  the 
mth  triangle  attached  to  node  j ,  and  the  column  is  defined  by  the  image  of 
the  £th  triangle  attached  to  node  i. 

Figure  4.3  illustrates  the  expansion  function  J^,  together  with  its  image 
J|.  Here  we  have  three  surfaces  attached  to  node  i.  It  is  seen  that  the  z 
component  of  the  expansion  function  is  the  same  on  each  surface  and  its  image. 
However  the  x  and  y  components  (which  are  parallel  to  the  ground  plane)  are 
reversed  in  direction.  Keeping  this  in  mind,  we  have  the  following  expres¬ 
sion  for  ZG (j , i) . 

t .  t . 

ZG ( j  ,  i)  =  -  f  Tjm  J.m  •  1  i  (K(jm,  i£)  -  K(jm,  i ’  Jl) ) 

(x(x  •  J_i?)  +  y(v  •  J^,))  +  (K(jm  ,  i£)  +  K(jm,  i'l))  z  (z  •  J^^)] 


But  from  (4-5)  we  have 


K(jm,  i '  i)  =  K(j'm,  iJ.)  . 


(4-17) 
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associated  with  node  i  and  its 


Hence, 
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ZG(j,i  ) 


T.  J . 
jm  -jm 


t  . 

1 

J  {  PG(jm,i£)  [x(x  •  J_i{)  +  ’  J^)] 

£  1 


+  [K(jm,  i£)  +  K(j'm,  i£)]  (z  •  j  )£} 

~ i£  — 


(4-18) 


Here  PG(jm,  i£)  is  an  element  of  matrix  PG  in  (4-9).  The  row  of  this 
element  is  defined  by  the  mth  triangular  surface  attached  to  node  j,  and 
the  column  is  defined  by  the  £th  triangle  attached  to  node  i. 

Hence,  we  use  the  electrostatic  matrix  PG  to  compute  the  magnetostatic 
matrix  ZG.  We  also  need  the  sum  K(jm,  i£)  +  K(j'm,  i£) ,  but  we  do  not  per¬ 
form  any  integration  to  find  this  sum,  because  both  of  its  terms  are  by¬ 
products  of  the  calculation  of  PG. 
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The  j th  element  of  VG  is  the  same  as  (2-49),  with  one  exception.  The 

exception  is  that  the  term  H^^r^^)  is  now  taken  to  mean  the  normal  component 

of  the  total  impressed  field  (i . e ., produced  by  the  true  source  plus  its 

image  )  at  r  ,  , 

-cjm 

Once  the  current  distribution  is  solved  for  using  (4-14),  the  induced 
magnetic  dipole  is  computed  using  (2-52).  This  dipole  replaces  the  body 
and,  strictly  speaking,  it  is  above  the  ground  plane.  Since  the  distance 
between  the  ground  plane  and  the  body  is  assumed  to  be  small,  we  place  the 
dipole  on  the  ground  plane  and  use  image  theory  to  remove  the  ground  plane. 

Figure  4.4  shows  the  induced  current  density  on  a  square  plate  which 
is  vertical  to,  and  1  cm  above,  the  ground  plane.  The  plate  is  1  meter 
by  1  meter  and  is  placed  in  the  yz  plane.  The  ground  plane  is  defined  by 
z  =  0.  The  incident  field  (in  the  presence  of  the  ground  plane  and  in  the 
absence  of  the  plate  )  is  a  z-polarized,  y-traveling,  plane  wave  (grazing 
incidence)  with  =  1  x  Ampere/meter.  Comparing  the  result  with  Fig.  3.15 

(  0  «*  90°  case  ) ,  which  represents  the  current  when  the  plate  is  in 
free  space,  we  observe  that  the  presence  of  the  ground  plane  destroys 


Fig.  4.4.  Current  on  a  square  plate  1  cm.  above  a  ground  plane.  No.  of  unknowns=25 
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the  symmetry  about  the  center  of  the  plate,  as  expected.  The  current  is 
bigger  near  the  edge  close  to  the  ground  plane.  The  computed  induced 
magnetic  dipole  is  M^nC*  =  -  0.455  x,  which  is  about  nine  percent  bigger  (in 
magnitude)  than  the  induced  dipole  moment  in  the  free  space  case  (-0.418  x) . 
Since  the  former  dipole  is  above  the  ground  plane,  we  can  conclude  the 
following  using  image  theory:  The  magnetic  (electric)  polarizability  of 
two  such  square  plates  (apertures  in  an  infinite  conducting  screen)  is 
larger  by  about  nine  percent  than  the  sum  of  their  individual  polariza¬ 
bilities  . 

Consider  a  small  cube  of  side  length  1  meter,  placed  a  distance  of 
0.5  meter  above  a  ground  plane.  The  ground  plane  is  at  z  =  0,  and  the 
top  face  of  the  cube  is  at  z  =  1.5  meters.  When  a  magnetic  dipole  is 
placed  tangentially  on  any  surface  of  the  cube,  we  use  the  pseudo-image 
method  described  in  Section  3.5.2  to  compute  the  induced  magnetic  dipole 
moment.  For  excitation  vector  elements,  we  use  (3-18).  However  in  this 
case  ElmP  of  (3-18)  is  produced  by  four  dipoles;  the  original  source  dipole, 
the  pseudo-image  of  the  original  dipole  just  inside  the  cube,  and  the  images 
of  these  two  dipoles  below  the  ground  plane.  Figure  4-5  illustrates  the 
problem  when  the  source  dipole  is  perpendicular  to  the  ground  plane.  When 
an  x-directed  magnetic  dipole  of  moment  =  1  is  placed  at  the  center  of 
the  top  face  of  the  cube  at  the  point  x  =  0,  y  =  (band  z  =  1.5,  the  com¬ 
puted  induced  dipole  is  given  by  M*nc^  =  0.18.  If  the  source  dipole  is 
placed  at  the  center  of  the  right  face  at  the  point  x  =  0,  z  =  l.O.and 
y  =  0.5,  the  induced  dipole  moment  is  M * n<^  =  0.2.  If  a  z-directed 
magnetic  dipole  of  moment  M  =  1  is  placed  at  the  center  of  the  right  face, 
the  induced  dipole  moment  is  M^nt*  =  0.21.  These  dipoles  represent  the  cube 


Fig.  4.5.  a)  Original  problem  :  A  magnetic  dipole  M  is  placed  tangentially 
on  the  right  face  of  a  cube  above  a  ground  plane,  b)  Equivalent  problem: 
is  the  pseudo- image  of  M.  M'  and  '  are  the  images  of  M  and  f 


respectively. 


above  the  ground  plane.  Since  a  z-directed  magnetic  dipole  has  an  image 
in  the  negative  z  direction,  we  conclude  that,  at  distances  large  compared 
to  the  size  of  the  cube,  the  scattered  field  will  be  smallest  for  the  last 
case  considered.  In  other  words,  to  prevent  radiation  from  sources 
inside  the  cube  to  the  outside  region,  a  narrow  slot  vertical  to  the  ground 
plane  is  better  than  a  slot  parallel  to  the  ground  plane. 

4.3  Conducting  Box  on  a  Ground  Plane 

Electronic  equipment  is  usually  contained  in  a  conducting  box  which  has 
several  apertures  for  input-output  connections  and  for  ventilation  purposes. 
It  is  desirable  to  keep  the  radiation  through  these  apertures  as  low  as 
possible.  To  investigate  the  effect  of  these  apertures,  test  measurements 
are  done  in  an  anechoic  chamber  with  the  box  usually  placed  on  a  ground 
plane.  In  this  section,  we  study  the  problem  of  a  box  on  a  ground  plane. 

The  box  is  assumed  to  have  a  narrow  slot,  and,  for  computational  purposes, 
this  slot  is  represented  by  a  magnetic  dipole  tangent  to  the  box.  Also, 
plane  wave  scattering  from  the  box  is  studied. 

4.3.1)  The  Electrostatic  Problem. 

The  electrostatic  problem  in  this  case  is  the  same  as  the  case  con¬ 
sidered  in  Section  4.1  with  two  exceptions.  Now  the  potential  is  zero, 
and  the  net  charge  is  unknown.  Hence  the  moment  equation  for  this  case  can 
be  written  as 
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where  all  quantities  appearing  in  this  equation  are  the  same  as  those 

in  (4-8). 
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Note  that  in  patching  the  box  surface  the  bottom  face  (the  part  of 
the  surface  on  the  ground  plane)  is  not  considered.  Hence,  the  patched 
surface  will  be  an  open  one.  For  a  rectangular  box,  the  boundary  of  this 
open  surface  is  a  rectangle  on  the  ground  plane.  The  image  of  the  patched 
surface  is  also  an  open  surface  with  the  same  boundary.  Hence,  the  patched 
surface  together  with  its  image  forms  a  closed  surface.  It  represents  a 
box  of  twice  the  height  of  the  original  box. 

4.3.2)  The  Magnetostatic  Problem 

The  moment  equation  for  this  problem  is  the  same  as  (4-14) ,  but  we 

have  the  following  modification  for  some  expansion  functions.  Since  the 

patched  surface  together  with  its  image  represents  a  closed  surface,  all  the 

nodes  are  considered  to  be  interior  nodes.  Hence,  we  associate  an  expansion 

function  with  each  node  but  one  .  The  image  of  a  node  on  the  boundary 

is  itself.  The  expansion  function  associated  with  such  a  node  cannot 

complete  a  loop  around  the  node.  However,  this  expansion  function  together 

with  its  image  forms  a  loop  around  the  node.  This  is  illustrated  in  Fig.  4.6. 

This  shows  the  expansion  function  associated  with  node  L,  which  is  assumed 

to  be  on  the  ground  plane.  An  expansion  function  associated  with  a  node  not 

on  the  ground  plane  is  the  same  as  the  one  illustrated  in  Fig.  2.2. 

Consider  a  cube  of  side  length  lm  placed  on  a  ground  plane  at  z  =  0. 

The  top  face  of  the  cube  is  at  z  =  lm,  and  the  right  face  is  at  y  =  0.5m. 

For  a  plane  wave  at  grazing  incidence  with  E*  =  377  (V/m)  and  Hx  =  1  (A/m),  the 

computed  induced  dipole  moments  are  =0.865*10  _z  and  11  =-1.32  21  •  (The 

scattered  fields  can  be  computed  by  placing  these  dipoles  on  the  ground  plane.) 

We  have  considered  three  magnetic  dipole  excitations  for  three 

different  boxes  on  the  ground  plane  .  First  ,  we  place  a  magnetic 

2 

dipole  with  dipole  moment  M  =  1  x  (  A-m  )  at  the  center  of  the  top 
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Fig.  4.6.  A  face  of  a  conducting  rectangular  box  together  with  its  image  in 
a  ground  plane. (The  ground  plane  is  at  z=0.)  The  solid  arrows  indicate  the 
expansion  function  associated  with  node  L.  The  dashed  arrows  represent  the 
image  of  the  expansion  function.  The  unprimed  capital  letters  denote  the 
nodes  and  the  primed  capital  letters  denote  the  images  of  the  nodes  in  the 
ground  plane. The  expansion  function  together  with  its  image  completes  a 
loop  around  the  node  L. 


face  of  the  box  .  This  is  referred  to  as  "first  excitation."  Secondly, 
we  place  M  =  1  x  at  the  center  of  the  right  face,  and  thirdly  we  place 
M  =  1  z  at  the  center  of  the  right  face.  The  boxes  considered  are 
i)  a  cube  (W  =  L  =  H  =  1  m) ,  ii)  a  tall  box  (W  =  L  =  0.9  m,  H  =  1.25  ir.) , 
and  iii)  a  short  box  (W  =  L  =  1.2m,  H=  0.7  m).  The  volumes  of  the 
three  boxes  are  almost  identical.  Table  4.1  summarizes  the  computed 
induced  magnetic  dipole  moments  obtained  for  these  three  boxes  with  each 
of  the  three  excitations .These  dipoles  radiate  in  the  presence  of  the  ground 
plane. 

Table  4.1.  The  induced  magnetic  dipole  moment  due  to  a  tangential 
magnetic  dipole  on  a  face  of  a  box. 


Short  box 


0.18  x 
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APPENDIX  A 

APPROXIMATION  OF  THE  LOW  FREQUENCY  SCATTERING  PROBLEM 


In  this  appendix  we  describe  how  to  approximate  the  problem  of  low 
frequency  electromagnetic  scattering  by  the  corresponding  electrostatic 
and  magnetostatic  problems.  The  procedure  is  explained  in  [1]  and  is 
summarized  here  for  quick  reference.  An  example  is  given  to  illustrate 
the  ideas  developed. 

Consider  a  perfectly  conducting  body  placed  in  an  unbounded  homogeneous 
medium  characterized  by  ( r, p) .  Let  the  body  be  illuminated  by  an  incident 
electromagnetic  field  (E1,^1)  whose  wavelength  is  much  larger  than  the 
maximum  linear  dimension  of  the  body.  The  problem  consists  of  determining 


s  „s. 


the  scattered  fields  (E  ,H  )  such  that 
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lim  r  x  (V  x  Ha) 


-  jkrH  =  0 


(A— 3) 

(A-4) 


where 

n  =  Ai/e  > 

k  =  U)  t/jjG  , 

0)  =  angular  frequency  of  the  incident  wave. 

S  denotes  the  surface  of  the  body,  and  n  is  the  unit  outward 
normal  vector  to  S. 

Note  that  in  (A-3)  the  boundary  condition  on  E'S  implies  the  boundary  con- 

s 

dition  on  H  . 


It  is  known  that  all  field  quantities  may  be  written  as 
convergent  power  series  in  k,  provided  k  is  sufficiently  small, 


that  is. 


'■*  CO 

E.S  * 1  =  [  (jk)m  ES  ’ 1  ;  J1S  ’ 1  =  l  (jk)m  HS,i 

m=0  “m  m=0  ~m 


(A 


The  mth  terms  in  these  expansions  may  be  found  by  substituting  (A-5) 
into  (A-l)-(A-3)  and  equating  the  coefficients  of  like  powers  of  k. 


We  find  that  the  first  terms,  ES  and  HS,  satisfy 

— o  — o  J 


V  x  E  =  0 


V  x  H  =  0 
—  — o 


ES  =  0 


HS  =  0 
— o 


n  x  t, 

—  — o 


=  -n  x  E 
—  — o 


n  •  H" 


=  -n  •  H 


In  addition  to  these  we  have  the  following  equations 


lim  r  E  =0 
— o 


lim  r  H  =0 
o 

f-KO 


(A- 6 


(A- 7) 


(A-8) 


(A- 9: 


These  equations  were  first  used  by  Rayleigh  [2]  and  subsequently  veri¬ 
fied  by  Stevenson  [20]  and  Kleinman  [21]. 

Note  that  the  boundary  conditions  in  (A-8)  are  now  independent.  Thi 
is  due  to  the  fact  that  equations  (A-6)  are  decoupled.  Equation  (A-6) 

g 

implies  that  there  exists  a  scalar  function  (J)  such  that 

o 

ES  =  -  V<j>S  (A-1C 

— o  —  o 

Then  the  left  equations  in  (A-7)  -  (A-9)  are  satisfied  if 
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2  s 

v  <r  =  o 

o 


(A-ll) 


-  n  x  V|S  =  -  n  x  e1 
~  ~°  S  ~  S 


(A-12) 


lim  r  —  =  0 

|*-X» 


(A-13) 


It  is  obvious  that  (A-11)-(A-13)  define  a  standard  electrostatic  Dirichlet 
potential  problem  for  the  surface  S  except  for  a  minor  detail.  The  excep¬ 
tion  being  that  the  boundary  condition  in  (A-12)  may  be  shown  to  specify 

s 

the  boundary  values  of  4>  within  an  arbitrary  additive  constant.  But  this 
constant  can  be  evaluated  by  employing  the  condition  that  there  is  no  net 
charge  on  the  surface  of  the  body,  i.e.. 


fl  •  E  ds  =  0 
—  — o 


(A-14) 


If  we  replace  (ES,  HS)  by  (E1,  H*)  in  (A-l)  and  use  (A-5),  we  have 


so  that 


V  x  E  =  0 
—  — o 


E1  =  -  V^)1 
— o  —  o 


where  d)1  is  a  scalar  function.  Then  (A-12)  implies  that 
o 


(A-15) 


(A-16) 


d)S  +  d)1  =  V  on  S 

o  o  o 


(A-17) 


where  V  is  a  constant, 
o 

We  can  think  of  the  problem  defined  by  (A-ll),  (A-13),  (A-14)  with 
s  s 

replaced  by  -V^q,  and  (A-17)  in  the  following  way.  A  perfectly  conducting 
isolated  body  (with  net  charge  Q  =  0)  is  placed  in  an  impressed  electrostatic 


potential  4>0  and  a  surface  charge  density  is  induced  on  the  surface 
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ci  the  body.  This  charge  produces  <j>®.  The  body,  being  a  perfect  conductor, 

is  an  equlpotential  surface  at  an  unknown  total  potential  V  ,  i.e., 

o 


4>s  (o  )  +  4)1  =  v 
0  0  o  o 


(A-18) 


This  represents  an  integral  equation  for  c  .  Once  it  is  solved  we  can 

o 

g 

obtain  <|>  at  an  arbitrary  point  r_  by 


o  (r’) 


i  ff  o  (.r  • ; 

*0  ■  ZF  I!  TFFF 
s’ 


(A-19) 


and  J:  is  given  by  (A-10) .  The  induced  electric  dipole  moment  can  be  ob¬ 
tained  by 


0  (r)r  ds 
o  -  — 


(A-20) 


In  a  similar  fashion  we  can  interpret  the  right  set  of  equations  in 

(A-6)  to  (A-9)  as  follows.  A  perfectly  diamagnetic  body  is  placed  in  an 

impressed  static  field  H^"  and  a  s  irface  current  J  is  induced  on  it, 

— o  — o 

g 

which  in  turn  produces  such  that  the  normal  component  of  the  total  mag¬ 
netic  field  is  zero  on  the  surface,  i.e.. 


n  •  (H^  +  h!(J  >)  =  0  on  S 
—  —0  —o  —o 


(A-21) 


This  represents  an  integral  equation  for  the  induced  surface  current  J  . 

o 

Additional  information  about  J  can  be  obtained  as  follows.  Consider  the 

— o 

continuity  equation 


„  Tmd  ,,  ind 

V  •  J  =  -  ike  o 
— s  — 


(A-22) 


where  c  =  1//JJE  is  the  speed  of  light  in  the  surrounding  medium,  and  J3 


and  O  are  respectively  the  current  and  charge  density  induced  on  the 
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surface  of  the  scatterer  due  to  the  incident  field  (E1,  H1)  V*  denotes  the 

—  —  *  — s 

surface  divergence  operator.  The  induced  sources  J^nd  and  a'*’n^  are  the 
sources  of  ES  and  HS  of  (A-l)  to  (A-4) .  Since  Jlnd  and  a^n<^  are  related 
(on  S)  to  the  total  fields  H1  +  jiS  and  E1  +  _ES,one  can  argue  that  J^11^  and 
alnd  can  be  expanded  into  convergent  series  of  the  form  given  in  (A-5)  if  k 
is  sufficiently  small.  Then  (A-22)  gives 


v  •  £(jk)»  T-nd  .  -  jkc  •)  (jk)"  oind 

m=0  m=0  m 

Equating  the  coefficients  of  like  powers  of  k,  we  obtain 


(A-23) 


V  •  Jlnd  =  0 
— s  — o 


(A-24) 


„  Tind  ind 

V  •  J,  =  -  C  0 
— s  — 1  o 


(A-25) 


Since  of  (A-21)  must  represent  .J  in  the  limit  as  k 


0  and  since 


Jind  =  lim  Jind 
k+0 


(A-26) 


we  conclude  that  J  =  Jdnd  and  hence 
— o  — o 


V  •  J  =  0 
— S  — o 


(A-27) 


This  is  the  property  of  stationary  currents.  That  is  why  we  call  the 
problem  defined  by  (A-21)  and  (A-27)  a  quasi-stationary  magnetic  problem 
or  simply  a  magnetostatic  problem.  Once  is  solved  for  using  (A-21)  and 

g 

(A-27)  the  scattered  field  H  can  be  computed  using  the  Biot-Savart  law. 

^o 

The  induced  magnetic  dipole  is  given  by 


Mlnd  =  \  ff  r  x  J  (r)ds 

2  JJ  ~  - 


(A-28) 
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Following  the  same  argument  that  led  us  to  (A-27)  we  conclude  that  Oq  of 

(A-18)  and  aoinc*  of  (A-25)  are  identical. 

At  distances  small  compared  with  the  wavelength,  the  scattered  field 

(E8 ,  Hs)  is  approximately  (E^,  H^)  where  _Eq  is  the  quasi-static  electric  field 

due  to  a  and  H  is  the  quasi-static  magnetic  field  due  to  J  .At  distances 

o—o  —0 

s  s 

that  are  appreciable  compared  with  the  wavelength,  (_E  ,  H  )  can  not  be 

approximated  by  (E  ,  H  )  because  (E  ,  H.)  is  not  a  propagating  field. 

— 0  O  -  O  “O 

s  s  ind 

Now,  (E  ,  H  )  is  the  field  radiated  by  J  .In  Appendix  B,  it  is  shown  that, 
at  distances  large  compared  with  the  maximum  dimension  of  a  small  surface,  the 
field  radiated  by  an  arbitrary  current  distribution  on  that  surface  is  the 
field  of  the  combination  of  the  electric  and  magnetic  dipole  moments  of  >J. 
Therefore,  at  distances  large  compared  with  the  maximum  dimension  of 
the  scatterer  ,  (Es,  HS)  is  the  field  radiated  by  the  combination  of  the 
electric  and  magnetic  dipole  moments  of  Jind  .  In  the  limit  as  k+0,  these 
dipole  moments  approach  pind  of  (A-20)  and  of  (A-28) . 

Notice  that  depending  on  the  structure  of  the  scatterer  and  the 
polarization  of  the  incident  wave  P*nd  and/or  M*nd  may  vanish.  For  example, 
if  a  plane  wave  impinges  on  a  plate  with  E1  (H1)  being  perpendicular  (tan¬ 
gential)  to  the  surface  then  P*nd  (M*nd)  is  zero. 

We  now  illustrate  the  ideas  developed  in  this  appendix  by  the  example 
of  plane  wave  scattering  from  a  small  conducting  sphere.  Figure  A.l  repre¬ 
sents  a  conducting  sphere  of  radius  a  illuminated  by  an  incident  plane  wave. 


Take  the  incident  wave  to  be  x-polarized  and  z-traveling,  that  is, 


i  -jkz  ^ 

E  =  V  x 


i  o  -jkz  -n 

=  -  e  J  v 


H  =  - 
-  n 


(A-29) 
(A- 30) 


where 
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n  =  /u  /e  =  120  7T  , 
o  o 


k  =  2~/\  =  wave  number  of  the  incident 


wave,  and  E  is  a  constant. 

o 


The  exact  results  are  taken  from  [4]. 

1 '  1 

jE  ,  00  sin  0  P  (cos  0)  P  (cos  0) 

__o  cos  <p  r  , _ n _ _  , _ n _ , 

0  n  ka  \  n  1  q (2) '  ,  3  .  A  r,(2)  „  ' 

n=l  H  (ka)  sm  0  H  (ka) 

n  n 

1  1 ' 

jE  .  ,  00  P  (cos  0)  sin  0  P  (cos  0) 

o  sin  <p  y  r  n  .  n 

4»  ‘  0  ka  A  3n  1  .  ~(2)'  H(2).,  , 

n=l  sin  0  H  (ka)  H  (ka) 

n  n 

jE  cos  00  /?\  i 

ES  - - — = —  T  b  n ( n+1 )  Hu;  (kr)  P  (cos  0) 

r  /t  .2  L  n  n  n 

(kr)  n=l 


(A-31) 


(A-32) 


(A-33) 


E  cos  (p  00  H^(kr)  P'*'(cos  0)  it 

Eq  - n -  I  [c  — - r — rji~ - jb  H'  (kr)P  (cos  0)sin  0] 

0  kr  4  1  n  sin  0  J  n  n  n  ‘ 

n=l 

(A-34) 

E  sin  4>  ™  ,,  (kr)P^(cos  0) 

E®  - -  l  [c  H'2  '  (kr)  sin  0  P  (cos  0)-jb  - •  o - ] 

d>  kr  L ,  n  n  n  n  sin  0 

n=l 

(A-35) 


a  =  - 


y  3  n ( 2n  +_1)  pl(cos  0) 
A  c.(2)\,  ,  VCOS 


n  2  '  ■'(2)' 

c(ka)  n=l  (ka) 

n 


(A-36) 


where 


J  and  J,  denote  the  0  and  (^-components  of  the  induced  surface 

0  <p 

current  density  respectively. 


E®,  E®,  E^  denote  the  r,  6. and  4>  components  of  the  scattered 
electric  field  respectively, 


O  is  the  induced  surface  charge  density,  and  c  is  the  speed  of  light. 


=  j  (2n+l) 
n  n(n+l) 


(A-37 


b  = 
n 


-  a 


J'  (ka) 
n 


nH(2)'(ka) 

n 


(A-38; 


J  (ka) 
n 


c  =  -  a 
n  n 

H(2)(k 

n 

J  (x)  =  X 

j  (x) 

n 

n 

xh(2) 

n 

(A-39; 


(A-40) 

(A-41) 


(2) 


Here»  Jn  and  hn  denote  the  spherical  Bessel  function  and  the  spherical 
Hankel  function  of  the  second  kind  respectively. 

The  distant  scattered  field  is  given  by 


0  kr-x”  kr 

JE 


jE  00  ,  P  (cos  6) 

—  e  J  cos  (J)  l  j  [b  sin  0  P  (cos  0)  -  c  5 — ]  (A-42) 

1  n  n  n  sxn  0  J 

n=l 


<p  kr^°°  kr 


-  ji,-  00  P  (cos  0) 

0  “Jkr  •  A  V  .nri.  n 

-  e  Sin  <t>  2.  J  [br 

n=l 


a' 


n  sin  0 


-  c^  sin  0  (cos  0)]  (A-43) 


Now  if  we  assume  that  k  is  small  such  that  (ka)  «  1,  then  we  have 


E1  =  x  Eq(1  -  jkz  +  (jk)2  ~  +  ...) 


.  E  2 

H  .  y  (1  -  jkz  +  (jkr  %-  +  ...) 


(A-44) 


Comparing  (A-44)  and  (A-5)  we  have 


E  =  x  E  =  -  V (-x  E  ) 


(A-45) 


From  (A-16)  and  (A-45)  we  have 


<£  =  -  x  E 

o  o 


Similarly  for  small  (ka)  we  obtain 


~  o  .  r  3  .  ,  7a  Q 

Jg  =  —  cos  d’f  -  —  +  (jk)  ~Y  cos  ft  +  . 


o  O  C 

J ,  =  —  sin  <t>[^r  cos  8  -  (jk)a  (77  +  7-  cos  26)  +  ...] 

<p  r i  z  Zb 


Li  r 

co  =  -2-  sin  0  cos  <f[3  -  (jk)  cos  0  +  ...]  I 

Comparing  (A-48)  -  (A-50)  with  the  series  representations  of  J*nd  and 
olnd  in  (A-23),  we  obtain 


J*nd  =  -r  — ~  f-  cos  <t>  8  +  cos  0  sin  (j)  $]  ,  (A- 

— 0  2  n  ~  — 

E  a 

jind  _  _o_  ^7  cog  ^  cqs  g  g  _  sin  (j,  (j-  +  -5.  cos  20)  j>] ,  (A- 

J.  n  3  /.  o 


ind  o  ~  >  a  1 

a  =  —  3  sin  0  cos  <p 
o  nc 


The  symbols  0  and  (p  denote  the  unit  vectors  in  the  0  and  (J)  directions 
respectively. 


Note  that 


Jind  =  0 


v  •  j: 

— S  —1 


in  agreement  with  (A-24)  and  (A-25) .  Observe  that  Oq  of  (A-53)  is  the 
static  charge  distribution  induced  on  the  sphere  due  to  the  impressed 
static  electric  field  given  in  (A-45).  Using  J^nd  and  o*nd  given  above 
we  compute  the  induced  dipole  moments  to  obtain 


field 


ind  ff  ind,  N 

P  =|  o  (r)  r 

-  jj  o  —  — 

sphere 


,  .  o  .  3 

ds  =  x  —  4tt  a 


so  that 


4tt  E  (ka) 

.  ind  .  .  o 
I*  =  JwP  =  x  j  - 5 - 

n  k 


(A-56) 


,ind  1 


,xnd.  .  , 
r  x  j  (r)ds 
—  — o  ~ 


sphere 


T  E 

~  3  o  ~ 

=  -  2tt  a  —  y 


so  that 


2^;  E 

K£  =  jojy  Mind  =  y  - (ka)^ 

°  jkZ 

Now  consider  the  distant  scattered  field.  For  small  (ka),(A-42)  and 
(A-43)  reduce  to 

-jkr  x 

Ee  Eo  (ka)  cos  ^  (cos  6  '  1/2) 


J  R.I.  o  , 

E*f  - — — — >E  -  (ka)  sin  4)  (—  cos  0 

t()  kr->°“  o  kr  T  2 


(A-57) 


(A-58) 


(A-59) 


The  field  given  by  (A-58)  and  (A-59)  is  the  far  field  of  the  combination 
of  the  dipoles  (A-56)  and  (A-57) .  This  is  in  agreement  with  the  fact 
that  the  scattered  far  field  of  a  small  body  can  be  expressed  in  terms 
of  an  electric  and/or  a  magnetic  dipole. 

Now  consider  the  expressions  given  in  (A-33)  -  (A-34) .  Assume 


ka  is  small  enough  such  that 


•j? 


96 


b 


n  ka-*-0 


n+1 _ r  2n  (n-I ) !  -I  2  (ka)2n+1 

n  Cn  ka-K)  (2n) !  n+1 


(A-60) 


Then  it  is  easy  to  show  that  for  any  kr  the  first  term  in  (A-33)  is  the 

g 

field  of  _It  in  (A-56) .  Hence  as  ka+0,  E^_  is  determined  by  _I?.  only. 

g 

Figure  A. 2  shows  the  scattered  field  Kc.  for  the  case  ka  =  0.2.  The  solid 
curves  represent  the  exact  solution  which  was  obtained  by  summing  the  first 
six  terms  of  (A-34) .  The  dashed  curves  represent  the  field  of  the  two  dipoles 
given  in  (A-56)  and  (A-57)  placed  at  the  origin.  The  fields  are  normalized 
with  respect  to  E^cos  t.  It  is  seen  from  the  figure  that  the  exact  fields  and 
the  dipole  fields  are  in  close  agreemei*.  even  at  a  distance  equal  to  twice 
the  radius.  We  believe  this  is  due  to  the  symmetrical  structure  of  the  sphere. 


A. 2.  Plane  wave  scattering  from  a  small  sphere. 


APPENDIX  B 


DIPOLE  REPRESENTATION  OF  THE  FIELDS  DUE  TO  ELECTRIC  AND  MAGNETIC 
CURRENTS  FLOWING  ON  A  SMALL  SURFACE 

Consider  a  time-harmonic  surface  density  of  electric  current  J  on 

the  surface  S  whose  maximum  dimension  D  is  appreciably  smaller  than  the 

wavelength.  This  electric  current  produces  an  electromagnetic  field 
s  s 

(El  » H  ) •  F°r  convenience,  the  coordinate  origin  is  chosen  such  that  the 

distance  r '  from  the  origin  to  any  point  on  S  is  less  than  D.  In  Appendix 

B,  it  is  shown  that,  whenever  the  distance  r  to  the  point  of  observation 

of  (E  ,H  )  is  large  compared  with  D,  (jl  ,H  )  reduces  to  the  electromagnetic 

field  of  the  combination  of  an  electric  dipole  and  a  magnetic  dipole,  both 

dipoles  being  located  at  the  origin.  From  inspection  of  this  dipole  repre- 

s  s 

sentation  of  the  electromagnetic  field  (E  ,H  )  due  to  the  electric  current  J, 
a  dipole  representation  of  the  electromagnetic  field  (ESm,HSm)  due  to  a 
magnetic  surface  current  .1™  on  S  is  easily  obtained. 

Jones  [22,  pp.  530-532]  has  shown  that  the  dipole  representation  is 
valid  under  the  more  restrictive  conditions  that  S  is  closed  and  perfectly 
conducting  and  that  r  is  not  only  large  compared  with  D  but  also  large 
compared  with  the  wavelength.  Not  assuming  that  S  is  perfectly  conducting 
but  assuming  that  S  is  closed,  Kleinman  [23,24],  in  a  long  but  straight¬ 
forward  presentation,  has  obtained  a  dipole  representation  for  the  electric 
field  ES  far  outside  S  due  to  an  electromagnetic  field  (E,H)  existing  on  S. 

g 

According  to  the  equivalence  principle  [4,  Sec.  3-5],  this  E  can  be  viewed 
as  the  electric  field  due  to  the  combination  of  the  electric  surface  current 
n  x  H  and  the  magnetic  surface  current  _E  v  n  where  £  is  the  unit  vector 
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which  points  outward  from  S.  In  order  to  obtain  his  final  result  [24, 
Eqs.  (16)— (18) ] ,  Kleinman  had  to  assume  that  r  is  large  compared  with  the 

wavelength . 

The  dipole  representations  of  Appendix  B  are  obtained  without 
assuming  that  S  is  perfectly  conducting  and  without  assuming  that  S  is 
closed.  These  representations  are  valid  whenever  r  is  large  compared 
with  D,  no  matter  how  r  compares  with  the  wavelength. 

In  search  of  the  dipole  representation  for  the  electric  field 
ES  due  to  the  electric  current  J,  we  write  [4,  Eqs.  (3-63)] 


_ES  =  t— r  (V(V  •  A)  +  k2A) 
J 


:(r ' ) 


A  i  i 

-  JJ 

s 


r-r 


ds ' 


(B-l) 


(B-2) 


In  (B-2),  _r'  is  the  radius  vector  (from  the  origin)  to  the  point  at 
which  the  differential  element  of  area  ds'  is  located  and  j:  is  the 
radius  vector  to  the  point  at  which  A  is  evaluated. 

Since  r'  is  on  S, 


r'  <  V 


(B-3) 


where  r'  is  the  length  of  the  vector  r_' .  If  it  is  assumed  that  r,  the 
length  of  satisfies 

r  »  D  (B-4) 


then 


r '  <<  r 


(B-3) 


Thanks  to  (B-5) ,  we  obtain 


I  r  -  r '  ■  -  r  -  r  •  r  ’ 


r  •  r 


- - T7  - :  -  (1  + 

r-r  r 


(B-6) 


1 


r 


( B— 7 ) 


vv 
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where  t_  is  the  unit  vector  in  the  £  direction.  Substitution  of  (B-6) 
and  (B-7)  into  (B-2)  gives 


"jkr 


e  '  ff 


-  ds' 


(B-8) 


Because  D  was  assumed  to  be  appreciably  smaller  than  the  wave¬ 
length,  it  is  evident  from  (B-3)  that 


kr'  «  1 


(B-9) 


Equation  (B-9)  authorizes  the  approximation 


ejW  ’  a'  :  1  t  jk  r  •  r' 


(E-10) 


Substituting  (B-10)  into  (B-8)  and  neglecting  the  second  order  term, 
we  obtain 

A  =  [}j  d(l')ds'  +  (jk  +  7>  jj  (£*!')  JCr’  >ds  ’  ]  (B-ll) 


The 


first  integral  in  (B-ll)  is  called  _I^. 


k  '  h 


J(r '  )ds ' 


(B— 1 2 ) 


The  dot  product  of  (B-12)  with  an  arbitrary  constant  vector  a  is 


—  *  — i  = 


ff 

a  •  J  (r ' )ds ' 

J  ~  ~  ~ 


( B— 1 3 ) 


Now, 


a  =  V ' (a  •  r') 


(B-14) 


where  V’  is  the  V  operator  with  respect  to  the  coordinates  of  r_' 
Substitution  of  (B-14)  into  (B-13)  gives 


S' 


a  •  I  =  | |  (7'(a  •  r'))  •  J(r')ds' 

-  -1  jj  -  - 


The  form  of  the  integrand  in  (B— 1 5 )  suggests  the  identity  [8,  Eq.  ( A— 1 ) ] 


W  •  V'<£ds'  =  - 


<J>(7'  *  W)ds'  +  (W  •  u,)dP' 
s  —  —  — b 


In  (B-16) ,  W  is  tangent  to  S,  and  •  is  the  surface  divergence  on  S 

-n.th  respect  to  the  coordinates  of  jr* .  If  S  is  closed,  the  integral  over 

C  is  to  be  deleted  from  (B-16).  If  S  is  open,  C  is  the  contour  that 

bounds  S,  d£'  is  the  differential  element  of  length  along  C,  and  u,  is 

— b 

the  unit  vector  tangent  to  S  and  normal  to  C.  The  direction  of  u,  is 

— b 

away  from  S.  Application  of  (B-16)  to  the  integral  on  the  right-hand 


side  of  (B-15)  gives 


f  f  f 

a  •  I,  =  -  |  (a  •  r')V'  •  J(r')ds'  +  I  (a  •  r')(J(r')  •  u^)d£'  (B-l 

—  —1  ii—  —  s -  - -  — b 


Since  no  line  charge  can  exist  on  C,  it  follows  that 


J(r')  •  u,  =0  onC  (B- 

-  — D 

In  view  of  (B-18)  and  the  fact  that  a  is  arbitrary,  (B-17)  implies  that 


h  =  - 


r'  V'  •  J(r ' )ds ’ 
—  s - 


The  second  integral  in  (B-ll)  is  called  I_^. 


l2  =  (r  •  r’)J(r’)ds' 


The  identity  formed  by  substituting  (B-12)  for  _I^  in  (B-19)  is  also 
valid  with  J(r')  replaced  by  (r  •  r')J(r').  As  a  result. 
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I-  =  -  [l  r'  V  *  ((r  •  r')J(r' 


)  J (r ' ) )ds ' 


(B-21) 


The  V^*  term  in  (E-21)  suggests  the  identity  [8,  Eq .  (A— 4) ] 


v'  •  (<t> w)  =  w  •  y’A  +  <tv*  •  w 

s  —  —  S  S  — 


( B— 2  2 ) 


In  (B-22),  W  is  tangent  to  S,  and  is  the  surface  gradient  on  S  with 
respect  to  the  coordinates  of  _r ' .  Application  of  (B-22)  to  the  integrand 
on  the  right-hand  side  of  (B-21)  gives 


I„  -  -  r'(J(r')  •  V* (r  •  r,))ds‘  -  r’(r  •  r’)  V’  •  J(r’)ds’  (B-23) 

—2  — -  s—  —  i  ,  ~  —  s - 


Since  J(r')  is  tangent  to  S,  the  surface  gradient  V'  in  (B-23)  can  be 

s 

replaced  by  the  ordinary  gradient  V'.  Using  (B-14)  with  a  replaced  by 
£  to  simplify  this  gradient,  we  obtain 


l2=~  r'(J(r')  •  r)ds '  -  r ’ (r  *  r')  •  J(r')ds' 


(B-24) 


Next,  ^2  expressed  as  half  of  the  sum  of  the  right-hand  sides  of 
(B-20)  and  (B-24). 


—2~\  ((£.  *  r_')J(r')  -  r'  (J(r')  •  r))ds'  -  jj  r '  (r  •  _r')V^  •  J(r')ds' 


(B-25) 


Thanks  to  the  vector  identity  [4,  last  of  Eqs.  (A— 1 3 ) ] ,  the  integrand 
of  the  first  integral  on  the  right-hand  side  of  (B-25)  is  -r  *  (_r*  x  .J) . 


As  a  result,  (B-25)  becomes 


-  2  —  '  i  J  1  *  J(_r  )ds '  ~  2  jJ  l'  (r  *  r')V^  •  J (r '  )ds* 
's  S 


The  equation  of  continuity  is 


Vs  *  J(£')  =  -ju>o(r’) 


where  a  is  the  surface  density  of  electric  charge  on  S.  In  view  of 
(B-27),  substitution  of  (B-19)  and  (B-26)  into  (B— 1 1 )  gives 


e  J  .  11 

A  =  [j«  I  £*  o( r')  (1  +  ±(jk  +  j)  (r 


r ' ) )ds '  - 


(jk  +  -)  r  x  I  r'  x  J  (r '  )ds '  ] 
r  —  l  i  — - 


It  is  evident  from  (B-5)  and  (B-9)  that 


|  (jk  +  4>(i  *  IT)  «  1 


As  a  result,  (B-28)  reduces  to 


..  1 


[jus  r'oCr^ds'  -  (jk  +  ^)r  x  ±  r'  x  J(r')ds'] 


If  kr^l,  (B-30)  reduces  to  the  electric  and  magnetic  dipole  terms  in 
(25,  Eq.  (14)  without  the  factor  Uo] .  Expression  (B-30)  can  be  rewritten 


ff  e"jkr  1  ff 

ja»  |J  r’  o(r')ds’  +  (V  x  -  Jj  r'  x  J(r’)ds’ 


Thanks  to  the  vector  identity  [4,  sixth  of  Eqs.  (A-14)],  (B-31)  becomes 
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where 


A  =  A  +  ~r—  V  x  f 

-  -  j'.'U 


A  = 


a-jkr 

4frr 


J“! 


ff 

J 

S 


r’  a (r 1 )ds ' 


-jkr 
e  J 

4:rr 


j03U 

2 


•f 

S 


r '  x  J (r ' )ds ' 


Substituting  (B-32)  for  A  in  (B-l),  we  obtain 


E®  =  -±-  (V(7  •  A)  +  k2A)  -  V  x  £ 

J'0£ 


(B-32) 


(B-33) 


(B-34) 


(B-35) 


According  to  [4,  first  of  Eqs.  (3-83)],  the  right-hand  side  of  (B-35) 
is  the  electric  field  of  the  combination  of  the  magnetic  vector  potential 
A  and  the  electric  vector  potential  _F.  Now,  (B-33)  can  be  expressed  as 


,  e 


-jkr 


-  “  —  4ur 


where 


I£  =  jw  [ 

J  J 


r’  G(r')ds' 


(B-36) 


(B-37) 


Since  A  is  a  magnetic  vector  potential,  l£  can  be  viewed  as  an  electric 
current  element  located  at  the  origin.  Equation  (B-34)  is  recast  as 


-jkr 

t  =  K£  % - 

—  —  4ttt 


(B-38) 


where 

K i  =  ff  r'  v  J(r '  )ds ' 

—  L  jj  - 

S 


(B-39) 


Since  £  is  an  electric  vector  potential,  K£  can  be  viewed  as  a  magnetic 
current  element  located  at  the  origin. 

Substituting  (B-36)  and  (B-38)  for  A  and  £  in  (B-35),  we  obtain 


105 


b 


ES  = 


l£  e 


-  j  kr 


(V(V‘ 


)  +  k 


2  11  e 


-  j  kr 


K£  e 


-  j  kr 


-  iioeVx- 


4tt  juie 

When  the  differential  operations  in  (B-40)  are  performed,  (B-40) 
expands  to 


-)  (B 


.s  e'jkr,.2 


rjk  .  JL 


-  =  4?i~~7  [k' XL  +  iff  +  -%)(3(l£-r)r-l£)  +  jwe(jk+^Xr  ^)] 

JU  r 

(B-41) 

g 

According  to  (B-35) ,  the  electric  field  E  due  to  J_  is  the  electric 
field  due  to  the  combination  of  A  and  £.  Hence,  the  magnetic  field  H  due 
to  £  must  be  the  magnetic  field  due  to  the  combination  of  A  and  £.  This 
magnetic  field  is  given  by  [4,  second  of  Eqs.  (3-83)] 


HS  =  t— -  (7(V  •  F)  +  k2F)  +  V  x  a 
j(4i  - 


Substituting  (B-36)  and  (B-38)  for  A  and  £  in  (B-42),  we  obtain 


(B-42) 


HS  =  t — t —  (V  (7 


K£  e 


-jkr 


K£  e'jkr  IP.  e"jkr 


-)  +  k 


2 


+  jLOpV  x 


-)  (B-43) 


Since  (B-43)  is  similar  to  (B-40),  it  is  evident  from  inspection  of  (B-41) 
that  (B-43)  expands  to 

”  j  kr 

Hs  =  -  [k2(?  xk£)  x  f  +  (JA  +  J-)(3(K£  •  r)r-K£)  -  jay (jk  +-)  (r  x  U)] 

—  4Trjcopr  —  —  ~  r  r  — - —  r  —  — 

(B-44) 

An  electric  charge  dipole  £  is  implicitly  defined  in  terms 
of  ££  by 

If  =  jwP  (B-45) 


-40) 


The  factor  jir  appears  in  (B-45)  because  electric  current  is  the  product 
of  jo)  with  electric  charge.  For  brevity,  £  is  called  an  electric  dipole. 


Alternatively,  £  can  be  called  the  electric  dipole  moment  of  the  electric 
current  £.  Comparison  of  (B-45)  with  (B-37)  gives 


P 


ft 

Jj 

S 


a(r')ds' 


(B-4l 


where  0  is  the  electric  charge  associated  with  J. 

A  magnetic  pole  charge  dipole  M  is  implicitly  defined  in  terms 
of  by 

Kf;  =  jwp  M  (B-4; 


Magnetic  pole  charge  is  magnetic  charge  divided  by  y.  The  factor  jwy 
appears  in  (B-47)  because  magnetic  current  is  the  product  of  jcoy  with 
magnetic  pole  charge.  For  brevity,  M  is  called  a  magnetic  dipole. 
Alternatively,  M  can  be  called  the  magnetic  dipole  moment  of  the  electric 
current  J.  Comparison  of  (B-47)  with  (B-39)  gives 


M  = 


[[  r'  x  j(r' 


2  j  J  L'  *  d(r')ds' 


(B-48 


Consider  a  surface  density  of  magnetic  current  £  on  S.  The 
electric  field  ESm  due  to  Jm  is  given  by  [4,  Eqs.  (3-4)  and  (3-5)] 


,sm 


=  -  V  x  F 


(B-49 


F 


r-r 


ds ' 


(B-50 


g 

On  the  other  hand,  the  magnetic  field  H  due  to  the  electric  current  J_ 
is  given  by  [4,  Eqs.  (3-4)  and  (3-5)] 

HS  =  7  x  a  (B-51 

where  A  is  given  by  (B-2) .  The  pair  of  equations  (B-49)  and  (B-50)  is 

stn 

similar  to  the  pair  of  equations  (B-51)  and  (B-2).  If  the  E^  in  (B-49) 


and  (B-50)  is  called  E  (J  )  and  if  the  H  in  (B— 5 1 )  and  (B-2)  is  called 


H  (J),  then 


ESm(jm)  =  -  HS(jm) 


(B-52) 


Directed  by  (B-52),  we  reverse  the  sign  of  the  right-hand  side  of  (B-44) 
and  replace  J  by  therein  to  obtain 

-jkr  i  ..  , 

Esm  =  JL. -  [k2(r  x  x;"1)  +  (^  +  4)(3(Um  •  r)r  -  I  lm)  +  jue  ( jk  +-)  (r  *  W m)  ] 

—  A'n'joJE'.r  —  —  —  r  — - —  r  —  — 

(B— 53 ) 


where  m  .  .,_/•/■ 

=  r’  x  Jm(r’)ds'  (B-54) 

S 

K'm  =  jojff  r'  om(r')ds’  (B-55) 

—  JJ  ~ 

S 

In  (B-55),  om  is  the  surface  density  of  magnetic  charge  associated 

.  ,  ,m 

with  J  . 


V’  •  Jm(r ' )  =  -  jwam(r ' ) 
s  — 


(B-56) 


Note  that  E^Sm  of  (B— 53)  is  Es  of  (B-41)  with  H,  and  K£  replaced 
by  If,1"  and  K£m,  respectively.  Hence,  the  magnetic  field  HSm  due  to  Jm 
can  be  obtained  from  inspection  of  (B-44). 


e-Jkr  r.  2 , 


,e.~  [kZ  (f  -  Kf”)  •  f  +  (Ji  +  i)O(K4"'-f)r-Kt"')-Juig0k+-r)(«xU“‘)] 
4TTjwvr  -  —  -  rr  — - —  r  “ 


(B-57) 


An  electric  charge  dipole  P™  is  implicitly  defined  in  terms  of 


I2m  by 


I 


Comparison  of  (B-58)  with  (B-54)  gives 


Pm  =  -  |  ||  r'  x  j"(r')ds' 


.m , 


,m 


A  magnetic  pole  charge  dipole  M  is  implicitly  defined  in  terms 


of  K2,m  by 


,,  ,m  .  m 
K*.  =  jwy  M 


Comparison  of  (B-60)  with  (B-55)  shows  that 


M 


m 


r’  am(r ' )ds ' 


(B 


-61) 
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